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of  the  McDonnell  Douglas  Corporation.  The  period  of  performance  was  May 
1987  to  April  1991.  Mr  Jerry  E.  Jenkins  was  the  Air  Force  project  engineer.  This 
*  report  supersedes  AFWAL-TR-86-3091  (Volume  II),  produced  under  Contract 

F33615-81-C-3617,  which  documents  Missile  Datcom  Rev  11/85,  12/88,  and 
7/89. 
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The  Aight  Dynamics  Directorate  is  committed  to  the  continuing 
development  of  Missile  Datcom.  This  development  is  dependent  to  a  large 
extent  on  user  feedback.  Questions  about  the  program  or  suggestions  for 
future  improvements  to  the  program  should  be  directed  to  Mr.  William 
Blake,  WL/FIGC,  Wright-Patterson  AFB,  Ohio  45433,  phone  (513)  255-6764. 
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BACKGROUND  TO  MISSILE  DATCOM 


Missile  Datcom  Rev  4/91  is  the  fifth  in  a  series  of  releases. 

In  the  late  1970's,  the  Tri-Service  Committee  on  Missiles  and  Projectile 
Aerodynamics  defined  the  need  for  a  Missile  Datcom  type  prediction  tool.  The  Air 
Force  was  chosen  as  the  lead  service  for  the  effort.  A  contract  was  subsequently 
awarded  to  the  McDonnell  Douglas  Astronautics  Company  (E33615-80-C-3605)  to 
recommend  specific  methods  for  inclusion  into  a  potential  computer  program  and 
identify  areas  where  further  work  was  needed.  The  final  report  from  this  effort, 
"Development  Feasibility  of  Missile  Datcom"  (AFWAL-TR-81-3130)  was  published  in 
October  1981. 

In  September  1981,  the  Missile  Datcom  Development  Contract,  (F33615-81-C-3617) 
was  awarded  to  the  McDonnell  Douglas  Astronautics  Company.  It  subdivided  the  effort 
into  four  distinct  phases.  The  initial  release  of  the  program  in  August  1984  represented 
the  "Phase  I"  interim  capability.  Cases  run  using  this  version  were  limited  to 
axisymmetric  bodies  with  no  more  than  eight  fins  total  (two  sets  with  up  to  four  fins 
each). 


The  second  release  of  the  program  (Rev  11/85)  represented  the  "Phase  IV" 
capability.  This  was  the  final  version  generated  under  Air  Force  Contract 

F33615-81  <1-3617.  It  added  capability  for  elliptic  bodies,  inlets  at  supersonic  speeds, 
dynamic  derivatives,  experimental  data  substitution,  and  coirfiguration  incrementing. 
It  also  increased  the  permissible  number  of  fins  to  thirty  two  (four  sets  with  up  to  eight 
fins  each).  Two  volumes  of  documentation  (User's  Manual  and  Program  Implementation 
Guide),  dated  November  1985,  were  printed. 

The  third  release  (Rev  12/88)  coincided  with  the  publishing  of 

AFWAL-TR-86-3091.  'Fhis  version  expanded  the  experimental  data  substitution  option 
and  dynamic  derivative  capability.  Errors  from  the  11/85  version  were  also  corrected. 
Volume  I  of  TR-86-3091  (Final  Report)  discusses  the  methods  selected  for  incorporation 
into  Missile  Datcom.  Volume  II  (User's  Manual)  is  an  updated  version  of  the  November 
1985  manual.  These  reports  are  available  from  DTIC  as  ADA-211086  and  ADA-210128 
respectively. 

The  fourth  release  (Rev  7/89)  added  little  new  capability,  its  primary  purpose  was 
to  correct  coding  errors,  expand  the  body-alone  dynamic  derivative  capability,  and 
modify  the  equivalent  angle-of-attack  formulation.  No  new  documentation  was 
published. 

Feedback  from  users  and  a  user  survey  revealed  many  areas  of  desired 
improvement  for  the  program.  Two  contracts  (F33615-86-C-3626  and  F33615-87-C-3604) 
were  awarded  to  Nielsen  Engineering  and  Research,  and  McDonnell  Douglas, 


respectively,  to  develop  new  or  improved  methodology.  The  4/91  version  of  the 
program  is  the  end  product  from  these  contracts.  The  major  changes  to  the  program  and 
associated  reference  material  are  listed  below; 

(a)  A  prediction  method  for  supersonic  plume/afterbody  interactions  has  been 
added  (refs.  1,2). 

(b)  The  method  for  inlet  aerodynamics  has  been  replaced  with  a  new  method 
(refs  1,3).  This  n^ethod  is  applicable  at  subsonic  and  supersonic  speeds, 
and  has  less  cumbersome  input  requirements. 

(c)  A  prediction  method  for  inlet  additive  drag  has  been  added  (refs  4,5). 

(d)  The  method  for  transonic  body  pressure  drag  (Mach  <  1.2)  has  been 
replaced  with  a  more  general  method,  applicable  to  any  nose  shape  (refs. 
1,6). 

(e)  A  prediction  method  for  body  protuberances  has  been  added. 

(0  The  method  for  fin  lateral  center  of  pressure  has  been  replaced  with  a  more 
general  method  (ref.  7,8). 

(g)  A  method  for  very  low  arpect  ratio  fins  has  been  added.  It  is  used  at  Mach 
numbers  greater  than  2.5  (ref.  9). 

(h)  A  prediction  method  for  high  aspect  ratio  fins  has  been  added  (ref.  10). 

(i)  The  fin  axial  fcwce  calculatior\s  have  been  revised. 

(]■)  The  partial  output  and  trim  output  have  been  expanded. 

(k)  The  airfoil  section  module  calculations  have  been  expanded. 

(l)  Coding  errors  have  been  corrected. 
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1.0  INTRODUCTION 


In  missile  preliminary  design  it  is  necessary  to  quickly  and 
economically  estimate  the  aerodynamics  of  a  wide  variety  of  missile 
configuration  designs.  Since  the  ultimate  shape  and  aerodynamic 
performance  are  so  dependent  upon  the  subsystems  utilized,  such  as  payload 
*  size,  propulsion  system  selection  and  launch  mechanism,  the  designer  must 

be  capable  of  predicting  a  wide  variety  of  configurations  accurately.  The 
^  fundamental  purpose  of  Missile  Datcom  is  to  provide  an  aerodynamic  design 

tool  which  has  the  predictive  accuracy  suitable  for  preliminary  design,  and 
the  capability  for  the  user  to  easily  substitute  methods  to  fit  specific 
applications. 
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2.0  PROGRAM  CAPABILITIES 


The  computer  code  is  capable  of  addressing  a  wide  variety  of 
conventional  missile  designs.  For  the  purposes  of  this  document,  a 
conventional  missile  is  one  which  is  comprised  of  the  following: 

•  An  axisymmetric  or  elliptically-shaped  body 

•  One  to  four  fin  sets  located  along  the  body  between  the  nose 
and  base.  Each  fin  set  can  be  comprised  of  one  to  eight 
identical  panels  attached  around  the  body  at  a  common 
longitudinal  position 

•  An  airbreathing  propulsion  system. 

To  minimize  the  quantity  of  input  data  required,  commonly  used 
values  for  many  inputs  are  assumed  as  defaults.  However,  all  program 
defaults  can  be  overridden  by  the  user  in  order  to  more  accurately  model  the 
configuration  of  interest. 

The  following  paragraphs  detail  the  configurations  that  can  be 
analyzed.  Later  paragraphs  describe  the  range  of  aerodynamic  coefficients  that 
can  be  predicted.  Finally,  the  program  constraints  are  discussed. 

2.1  Addressable  Configurations 

The  following  configurations  can  be  analyzed: 

•  Circular  or  elliptically-shaped  cross  section  bodies,  with  or 
without  airbreathing  inlets 

•  Fin  alone  (1  to  8  panels  attached  at  the  root) 

•  Body  and  up  to  4  fin  sets  (1  to  8  panels  in  each  fin  set) 

•  The  body  and  fin  set  configurations  with  deflected  fins 

Certain  restrictions  exist  due  to  method  limitations  and  are 
summarized  in  the  following  paragraphs. 

2.1.1  Axisymmetric  or  Elliptical  Bodies  -  Methodology  is  incorporated  that 
permits  analysis  of  the  configuration  components  summarized  in  Table  1. 
Due  to  the  types  of  methods  selected  restrictions  also  apply  to  the  manner  in 
which  these  components  are  joined  to  form  a  complete  configuration: 
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•  Subsonic/transonic  speeds  -  The  aerodynamic  methods 
assume  that  the  body  is,  as  a  minimum,  composed  of  a  nose- 
cylinder  combination.  The  afterbody  (boattail  or  flare)  is 
optional,  but  if  used,  it  must  be  attached  to  a  cylindrical  center 
body  whose  length  is  at  least  four  body  diameters;  this 
restriction  minimizes  nose  flow  field  coupling  over  the 
afterbody.  If  an  afterbody  is  specified  it  must  not  be 
cylindrical,  e.g.,  the  base  diameter  must  be  different  than  the 
centerbody  diameter.  Table  2  summarizes  the  other 
restrictions  on  the  configurations. 

•  Supersonic  speeds  -  The  aerodynamic  methods  used  are  not 
restricted  to  nose-cylinder  combinations.  Any  arbitrary  radii 
distribution  can  be  defined  since  theoretical  techniques  are 
employed  at  Mach  numbers  above  1.2.  Care  should  be  taken 
to  avoid  introducing  unexpected  corners  into  the  contour.  If 
the  contour  has  any  concaved  regions  the  marching  may  fail 
due  to  shock  impingement  on  the  body  as  it  starts  to  curve 
out. 

2.1.2  Panels  -  The  prograim  will  accept  inputs  to  describe  most  airfoil  sections 
or  planforms.  Certain  assumptions  and  limitations  are  made  and 
summarized  in  the  following  paragraphs. 

2. 1.2.1  Airfoil  Section  -  The  program  will  accept  virtually  any  symmetrical 
airfoil  section  or  NASA  subsonic  cambered  section.  The  airfoil  section  can  be 
defined  using  a  NACA  designation  or  by  supplying  the  coordinates  of  the 
section.  Circular  arc,  hexagonal,  or  diamond  shaped  sections  can  also  be 
specified,  A  symmetric  hexagonal  cross-section  is  the  default;  its  shape  is 
computed  using  the  planform  inputs.  Hence,  explicit  definition  of  the  airfoil 
section  is  optional.  Although  cambered  airfoil  sections  can  be  input,  their  use 
in  the  code  is  currently  limited  to  subsonic  applications. 

2.1  2.2  Planform  -  Each  set  of  fins  may  be  comprised  of  up  to  eight  separate 
panels.  It  is  assumed  that  each  panel  is  geometrically  identical.  Although 
planforms  may  be  described  by  up  to  10  separate  pieces  or  sections,  an 
equivalent  straight-tap>ered  panel  is  computed  and  used  at  all  speeds.  There  is 
no  capability  to  specify  a  panel  with  outboard  dihedral. 

2.1.3  Airbreathing  Inlets  -  Both  axisymmetric  and  two-dimensional 
airbreathing  inlet/diverter  combinations  can  be  defined.  Up  to  20  identical 
inlets  can  be  positioned  around  the  body  at  arbitrary  angles.  Vehicles  with 
inlets  can  be  analyzed  at  all  speeds. 

2.2  Types  of  Data  Computed 
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2.2.1  Aerodynamics  -  The  program  computes  the  following  aerodynamic 
parameters  as  a  function  of  angle  of  attack  for  each  configuration: 

Normal  Force  Coefficient 
Lift  Coefficient 
Pitching  Moment  Coefficient 

Center  of  Pressure  in  calibers  from  the  moment 
reference  center 
Axial  Force  Coefficient 
Drag  Coefficient 
Side  Force  Coefficient 
Yawing  Moment  Coefficient 
Rolling  Moment  Coefficient 

Normal  force  coefficient  derivative  with  angle  of 
attack 

Pitching  moment  coefficient  derivative  with  angle 
of  attack 

Side  force  coefficient  derivative  with  sideslip  angle 
Yawing  moment  coefficient  derivative  with 
sideslip  angle 

Rolling  moment  coefficient  derivative  with 
sideslip  angle 

The  derivative  output  can  be  in  degrees  or  radians.  Partial  output 
results,  which  detail  the  components  used  in  the  calculations,  are  also 
optionally  available. 

It  should  be  noted  that  the  drag  force  (and  drag  coefficient)  is  different 
between  the  wind  and  stability  axes  systems  if  the  missile  body  is  at  a  sideslip 

angle  (P)  to  the  wind.  However,  wind  axis  drag  and  stability  axis  drag  are  the 
same  at  zero  sideslip.  In  Missile  Datcom,  drag  force  methods  are  assumed  to 
be  in  the  stability  axes  system  and  axial  force  methods  are  assumed  to  be  in 
the  body  axes  system  unless  otherwise  noted. 

The  program  has  the  capability  to  perform  a  static  trim  of  the 
configuration,  using  any  fin  set  for  control  with  fixed  incidence  on  the  other 
sets.  The  two  types  of  aerodynamic  output  available  from  the  trim  option  are 
as  follows: 

•  Unirimmed  data  -  Each  of  the  aerodynamic  force  and 
moment  coefficients  are  printed  in  a  matrix,  which  is  a 
function  of  angle  of  attack  and  panel  deflection  angle.  This 
output  is  optional. 


Cn 

Cl 

Cm 

Xcp 

Ca 

Co 

Cy 

Cn 

C| 

Cno 

Cma 


Cyp 

Cnp 

Cip 
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•  Trimmed  data  -  The  trimmed  aerodynamic  coefficients,  and 
trim  deflection  angle,  are  output  as  a  function  of  angle  of 
attack. 

2.2.2  Geometry  -  All  components  of  the  configuration  have  their  physical 
properties  calculated  and  output  for  reference  if  requested.  All  data  is  supplied 
in  the  user  selected  system  of  units. 

2.2.3  Other  -  The  reference  area  and  reference  length  are  user  defined.  The 
user  may  optionally  select  to  print  the  calculated  body  or  fin  pressure 
coefficient  distributions  at  supersonic  speeds.  Outputs  of  the  partial 
aerodynamic  results  and  a  summary  of  method  extrapolations  are  also 
optionally  available. 

2.3  Operational  Considerations 

The  code  has  been  written  to  conform  to  the  coding  standards  for  the 
American  National  Standards  Institute  (ANSI)  Standard  FORTRAN  X3.9- 
1978,  often  referred  to  as  FORTRAN  V  or  FORTRAN  77. 

There  are  only  four  exceptions  to  the  ANSI  standards  used  in  the 
computer  code: 

•  Transfer  on  end-of-file  -  The  FORTRAN  IV  statement 
IF(EOF(L1NIT).NE.O)  is  used  in  the  standard  CDC  compatible 
code.  The  FORTRAN  V  statement  END=label  in  the 
FORTRAN  READ  statements  is  incorporated  in  the  code,  but 
is  inactive  in  the  CDC  code. 

•  NAMELIST  -  The  use  of  nam.elist  for  input  and  output  (I/O) 
is  used.  Although  this  appears  to  be  a  violation  of  FORTRAN 
IV,  it  is  really  not  since  a  namelist  emulator  has  been  written 
for  the  Missile  Datcom  code  using  FORTRAN  IV. 

•  Mathematical  functions  -  A  few  mathematical  functions  are 
not  considered  "standard,"  such  as  the  trigonometric  tangent. 
Standard  FORTRAN  equivalents  for  these  functions  are 
available  on  request. 

•  PROGRAM  card  -  The  code  was  developed  on  the  Control 
Data  Corporation  CYBER  computers.  This  system  requires 
that  the  first  card  of  the  main  routine  be  a  PROGRAM 
statement.  An  IBM  or  VAX  compatible  version  of  the  code  is 
also  maintained  which  has  a  different  format  for  the  program 
card. 
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3.0  INPUT  DEFINITION 


Inputs  to  the  program  are  grouped  by  "case”.  A  "case"  consists  of  a  set 
of  input  cards  which  define  the  flight  conditions  and  geometry  to  be  run. 
Provisions  are  made  to  allow  multiple  cases  to  be  run.  The  successive  cases 
can  either  incorporate  the  data  of  the  previous  case  (using  the  input  card 
SAVE)  or  be  a  completely  new  configuration  design.  The  SAVE  feature,  for 
example,  permits  the  user  to  define  a  body  and  wing  (or  canard)  configuration 
in  the  first  case  and  vary  the  tail  design  for  subsequent  cases. 

The  scheme  used  to  input  data  to  the  computer  program  is  a  mixture  of 
namelist  and  control  cards.  This  combination  permits  the  following; 

•  Inputs  are  column  independent  and  can  be  input  in  any 
order. 

•  All  numeric  inputs  are  related  to  mnemonic  (variable) 
names. 

•  Program  input  "flags"  are  greatly  reduced.  Required  "flags" 
are  identified  by  a  unique  alphabetic  name  which 
corresponds  to  the  option  selected. 

The  program  includes  an  error  checking  routine  which  scans  all  inputs 
and  identifies  all  errors.  This  process  is  a  single-pass  error  checking  routine; 
all  errors  are  identified  in  a  single  "run".  In  addition,  the  program  checks  for 
necessary  valid  inputs,  such  as  a  non-zero  Reynolds  number.  In  some  cases, 
the  code  will  take  corrective  action.  The  type  of  corrective  action  taken  is 
summarized  later  in  this  section. 

Flexibility  has  been  maintained  for  all  user  inputs  and  outputs.  The 
following  summarize  the  program  generality  available: 

•  The  units  system  can  be  feet,  inches,  meters  or  centimeters. 

The  default  is  feet. 

•  Derivatives  can  be  expressed  in  degree  or  radian  measure. 
Degree  measure  is  the  default. 

•  The  body  geometry  can  be  defined  either  by  shape  typo  or  by 
surface  coordinates. 

•  The  airfoil  can  be  user  defined,  NACA,  or  supersonic  shaped 
sections.  The  NACA  sections  are  defined  using  the  NACA 
designation.  A  hexagonally  shaped  supersonic  section  is  the 
default. 
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•  The  configuration  can  be  run  at  a  fixed  sideslip  angle  and 
varying  body  angle  of  attack,  or  a  fixed  aerodynamic  roll  angle 
and  varying  total  angle  of  attack. 

•  The  flight  conditions  can  be  user  defined,  or  set  using  a 
Standard  Atmosphere  model.  The  capability  to  define  wind 
tunnel  test  conditions  as  the  flight  conditions  is  also 

available.  The  default  flight  condition  is  zero  altitude.  * 

3.1  Namelist  Inputs 

The  required  program  inputs  use  FORTRAN  namelists.  Missile 
Datcom  is  similar  to  other  codes  which  use  the  namelist  input  technique,  but 
differ  as  follows; 

•  Namelist  inputs  are  column  independent,  and  can  begin  in 
any  column  including  the  first.  If  a  namelist  is  continued  to  a 
second  card,  the  continued  card  must  leave  column  1  blank. 

Also,  the  card  before  the  continued  card  must  end  with  a 
comma.  The  last  usable  column  is  number  79  if  column  1  is 
used,  and  column  80  if  column  1  is  blank. 

•  The  same  namelist  can  be  input  multiple  times  for  the  same 
input  case.  The  total  number  of  namelists  read,  including 
repeat  occurrences  of  the  same  namelist  name,  must  not 
exceed  300. 

The  three  namelist  inputs 

$REFQ  SREF=1.,$ 

$REFQ  LREF=2.,$ 

$REFQ  ROUGH=0.001,$ 

are  equivalent  to 

$REFQ  SREF=1.,LREF=2.,ROUGH=0.001,$ 

•  The  last  occurrence  of  a  namelist  variable  in  a  case  is  the 
value  used  for  the  calculations. 

The  three  namelist  inputs 

$REFQ  SREF=1.,$ 

SFLTCON  NMACH-2.,MACH=1 .0,2.0,$ 

$REFQ  SREF=2.,$ 
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are  equivalent  to 


t 


$REFQ  SREF=2.,$ 

SFLTCON  NMACH=2.,  MACH=1.0,  2.0,$ 

•  The  namelists  can  be  input  in  any  order. 

•  Only  those  namelists  required  to  execute  the  case  need  be 
entered. 

•  Certain  hollerith  constants  are  permitted.  They  are 
summarized  in  Table  3.  Note  that  any  variable  can  be 
initialized  by  using  the  constant  UNUSED;  for  example, 
LREF=UNUSED  sets  the  reference  length  to  its  initialized 
value. 


All  Missile  Datcom  namelist  inputs  are  either  real  numbers  or  logical 
constants.  Integer  constants  will  produce  a  nonfatal  error  message  from  the 
error  checking  routine  and  should  be  avoided. 

The  namelist  names  have  been  selected  to  be  mnemonically  related  to 
their  physical  meaning.  The  ten  namelists  available  are  as  follows: 


Namelist 

SFLTCON 

SREFQ 

SAXIBOD 

SELLBOD 

SPROTUB 

SFINSETn 

SDEFLCT 

$TRIM 

SINLET 

SEXPR 


Inputs 

Flight  Conditions  (Angles  of  attack, 
Mach  numbers,  etc.) 

Reference  quantities  (Reference  area, 
length,  etc.) 

Axisymmetric  body  definition 
Elliptical  body  definition 
Protuberance  information  and 
geometry 

Fin  descriptions  by  fin  set  (n  is  the  fin 
set  number;  1,  2, 3  or  4) 

Panel  incidence  (deflection)  values 
Trimming  information 
Inlet  geometry 
Experimental  data 


Each  component  of  the  configuration  requires  a  separate  namelist 
input.  Hence,  an  input  case  of  a  body-wing-tail  configuration  requires  at  least 
one  of  each  of  the  following  namelist  inputs,  since  not  all  variables  have 
default  values  assigned: 


SFLTCON  to  define  the  flight  conditions 
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$AXIBOD  or  $ELLBOD  to  define  the  body 
$FINSET1  to  define  the  n\ost  forward  fin  set 

$FINSET2  to  define  the  first  following  fin  set 

$FINSET3  to  define  the  second  following  fin  set 

$FINSET4  to  define  the  third  following  fin  set 

The  following  namelists  are  optional  since  defaults  exist  for  all  inputs: 

$REFQ  to  define  the  reference  quantities 

$PROTUB  to  define  protuberance  option  inputs 

SDEFLCT  to  define  the  panel  incidence 

(deflection  angles) 

$TRIM  to  define  a  trim  case 

$INLET  to  define  inlet  geometry 

$EXPR  to  define  experimental  input  data 

Defaults  for  all  namelists  should  be  checked  to  verify  the  configuration 
being  modeled  does  not  include  an  unexpected  characteristic  introduced  by  a 
default. 

The  following  sections  describe  each  of  the  namelist  inputs.  Each 
section  is  accompanied  by  a  figure  which  summarizes  the  input  variables, 
their  definitions,  and  units.  Since  the  system  of  units  can  be  optionally 
selected,  the  column  "Units"  specifies  the  generic  system  of  units  as  follows: 

L  Units  of  length;  feet,  inches,  centimeters  or  meters 
F  Units  of  force;  pounds  or  Newtons 
deg  Units  of  degrees;  if  angular,  in  angular  degrees;  if 
temperature,  either  degrees  Rankine  or  degrees 
Kelvin 

sec  Units  of  time  in  seconds 

Expo  ients  are  added  to  modify  the  above.  For  example,  means  units 
of  length  squared,  or  area.  Combinations  of  the  above  are  ^so  used  to  specify 
other  units.  For  example,  F/L^  means  force  divided  by  area,  which  is  a 
pressure. 

Since  it  is  difficult  to  discern  the  difference  between  the  number  zero 
"0"  and  the  alphabetic  letter  "O",  it  should  be  noted  that  none  of  the  namelist 
or  namelist  variable  names  contain  the  number  zero  in  them.  In  general,  the 
number  zero  and  the  letter  "O"  are  not  interchangeable  unless  so  stated. 

The  program  ascertains  the  configuration  being  modeled  by  the 
presence  of  each  component  namelist,  even  if  no  data  is  entered.  The 
following  rules  for  namelist  input  apply: 
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•  Do  not  include  ^  unless  it  js  required.  Once  read,  the 

presence  of  a  n^rfieiist  (and/  hence,  a  configuration 
component)  can  oh^y  ^  removed  using  the  DELETE  control 
card  in  a  subseqnet'l  ca^e.  Simply  sotting  all  variables  to  their 
initialized  valUe^  Will  not  rentove  the  configuration 
component. 

•  Do  not  include  a  Variable  within  a  namelist  unless  it  is 
required.  Prograrh  adions  are  often  determined  from  the 
number  and  t)^s  input  providecl- 

•  Do  not  over-sp^cii^y  ihe  geometry-  User  inputs  will  take 

precedence  ovef  calculations.  Inputs  that  define  a 

shape  that  is  physi<^ally  impossibly  will  be  used  as  specified. 

The  program  doys  Vot  ’'fix-up"  inconsistent  or  contradictory 
inputs, 

3.1.1  NAMELIST  FLTCQM,  -  hight  Conditions 

This  namelist  defines  flight  conditions  to  be  run  for  the  case.  The 
program  is  limited  to  no  th^ui  20  angles  c>f  attack  and  20  Mach  numbers 
per  case  at  a  fixed  sideslip  aAgle.  aerodynanu^  roll  angle,  altitude,  and  panel 
deflection  angle.  Therefofe/  a  "trase"  is  defined  as  a  fixed  geometry  with 
variable  Mach  number  and  of  attack- 

The  inputs  are  giv^n  l^igure  1-  Tfiefe  are  two  ways  in  which  the 
aerodynamic  pitch  and  yaw  ahgle^  can  be  dcfiped: 

•  Input  ALPHA  ahd  If  BETA  is  input  and  PHI  is  not,  it  is 
assumed  that  tho  b^dy  a^is  angles  of  attack  (a)  and  sideslip 
angles  (p)  are  dcfin^. 

•  Input  alpha  and  pHl.  If  PHI  is  input  and  non-zero,  it  is 
assumed  that  ALPI^a  the  total  angle  of  attack  (a)  and  PHI 
is  the  aerodynarhic  •'nfi  ^gle  (<!’). 

•  Input  ALPHA,  bET^  and  hHI.  The  value  for  BETA  is  ignored 
if  PHI  is  non-zefO. 

As  a  minimum  the  folinw'ng  variables  must  be  defined: 

NALPHA  Auhlpe^.  of  angles  of  attack  to  run 

(MAlFHA  2  2) 

ALPHA  angle  of  attacl;  schedule  (matching 

NAlpftA) 
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NMACH 


number  of  Mach  numbers  or  speeds 
(NMACH  ^  1) 

MACH  or  VINF  Mach  number  or  speed  schedule  (matching 
NMACH) 

The  REN,  TEMF  and  PESIF  data  must  correspond  to  the  MACH  or  VINF 
inputs.  The  ALPHA  and  MACH  dependent  data  can  be  input  in  any  order;  the 
code  will  sort  the  data  into  ascending  order. 

Reynolds  number  is  always  required.  Three  types  of  inputs  are 
permitted  to  satisfy  the  Reynolds  number  requirement: 

•  Specify  Reynolds  number  per  unit  length  using  REN 

•  Specify  the  altitude  using  ALT,  and  the  speed  using  MACH  or 
VINF  (Reynolds  number  is  computed  using  the  Standard 
Atmosphere  model) 

•  Specify  pressure  and  temperature  using  PINF  and  TINF,  and 
the  speed  using  MACH  or  VINF  (typical  of  data  available 
from  a  wind  tunnel  test) 

User  supplied  data  will  take  precedence  over  program  calculations. 
Hence,  the  user  can  override  any  default  or  Standard  Atmosphere  calculation. 
The  default  condition  is  sea-level  altitude  (ALT=0.)  if  the  wrong  combination 
of  inputs  are  provided  and  the  Reynolds  number  cannot  be  calculated. 

3.1.2  NAMELIST  REFQ  -  Reference  Quantities 

Inputs  for  this  namelist  are  optional  and  are  defined  in  Figure  2.  A 
vehicle  scale  factor  (SCALE)  permits  the  user  to  input  a  geometry  that  is 
scaled  to  the  size  desired.  This  scale  factor  is  used  as  a  multiplier  to  the  user 
defined  geometry  inputs;  it  is  not  applied  to  the  user  input  reference 
quantities  (SREF,  LREF,  LATREF).  If  no  reference  quantities  are  input,  they 
are  computed  based  upon  the  scaled  geometry.  XCG  is  input  relative  to  the 
origin  (X=0)  and  is  scaded  using  SCALE. 

In  lieu  of  specifying  the  surface  roughness  height  ROUGH,  the  surface 
Roughness  Height  Rating  (RHR)  can  be  specified.  The  RHR  represents  the 
arithmetic  average  roughness  height  variation  in  millionths  of  an  inch. 
Typical  values  of  ROUGH  and  RHR  are  given  in  Table  4. 

3.1.3  NAMELIST  AXIBOD  -  Axisymmetric  Body  Geometry 

An  axisymmetric  body  is  defined  using  this  namelist.  The  namelist 
input  variables  are  given  in  Figures  3a  and  3b  and  a  sketch  of  the  geometric 
inputs  are  given  in  Figure  4.  The  body  can  be  specified  in  one  of  two  ways: 
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OPTION  1:  The  geometry  is  divided  into  nose,  centerbody,  and 
aft  body  sections.  The  shape,  overall  length,  and  base  diameter 
for  each  section  are  specified.  Note  that  not  all  three  body 
sections  need  to  exist  on  a  configuration;  for  example,  a  nose- 
cylinder  configuration  does  not  require  definition  of  an  aft  body. 

OPTION  2:  The  longitudinal  stations  and  corresponding  body 
radii  are  defined,  from  nose  to  tail.  This  option  should  only  be 
selected  if  the  Mach  number  is  greater  than  1.2. 


The  program  uses  the  input  value  for  NX  to  determine  which  option  is 
being  used.  If  NX  is  not  input  then  Option  1  inputs  are  assumed.  If  both 
shapes  and  body  coordinates  (Options  1  and  2)  are  used,  the  body  coordinate 
information  will  take  precedence.  NX  can  be  set  to  its  initialized  value  (to 
simulate  the  variable  as  not  input)  by  specifying  '’NX=UNU5ED". 


It  is  highly  recommended  that  Option  1  be  used  when  possible.  The 
program  automatically  calculates  the  body  contour  based  upon  the  segment 
shapes  using  geometry  generators.  Hence,  more  accurate  calculations  are 
possible.  Even  when  Option  2  is  used,  appropriate  Option  1  inputs  should  be 
included.  This  identifies  where  the  code  should  insert  break  points  in  the 
contour.  If  these  parameters  are  not  input,  they  are  selected  as  follows: 


LNOSE 

DNOSE 

LCENTR 

DCENTR 

LAFT 

DAFT 

DEXIT 


Length  of  the  body  segment  to  where  the 
radius  first  reaches  a  maximum 
The  diameter  at  the  first  radius  maximum 
Length  of  the  body  segment  where  the  radius 
is  constant 

Diameter  of  the  constant  radius  segment 
The  remaining  body  length 
Diameter  at  the  base 

Not  defined  (implies  that  base  drag  is  not  to 
be  included  in  the  axial  force  calculations) 


If  DEXIT  is  not  input,  or  set  to  UNUSED,  the  base  drag  computed  for 
the  body  geometry  will  not  be  included  in  the  final  computed  axial  force 
calculations.  To  include  a  "full"  base  drag  increment,  a  zero  exit  diameter 
must  be  specified  (DEXIT=0.). 

If  body  coordinates  are  input  using  the  variables  NX,  X,  R,  and 
DISCON,  and  the  nose  is  spherically  blunted,  the  geometry  must  be 
additionally  defined  using  the  following: 


•  BNOSE  must  be  specified  (even  if  zero) 


•  TRUNC  must  be  set  to  .FALSE. 
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•  The  first  five  (5)  points  in  the  X  and  R  arrays  must  lie  on  the 
spherical  nose  cap  [i.e.,  X(l),  X(2),  X(3),  X(4),  X(5),  Rvl),  R(2), 

R(3),  R(4),  and  R(5)  are  spherical  cap  coordinates] 

The  following  summarizes  the  input  generality  available: 

•  X(l)  does  not  have  to  be  0.0;  an  arbitrary  origin  can  be 
selected. 

•  Five  shapes  can  be  specified  by  name; 

CONICAL  (CONE)  -  cone  or  cone  frustrum  (default  for 
boattails  and  flares) 

OGIVE  -  tangent  ogive  (default  for  noses) 

POWER  -  power  law* 

HAACK  -  L-V  Haack  (length-volume  constrained)* 

KARMAN  -  von  Karman  (L-D  Haack;  length-diameter 
constrained)* 

•  If  DAFT<DCENTR  the  afterbody  is  a  boattail.** 

•  If  DAFr>DCENTR  the  afterbody  is  a  flare.** 

•  If  LAFT  is  not  input,  aft  body  (boattail  or  flare)  does  not  exist. 

•  applies  to  noses  only 

•  *  DAFT  must  not  be  equal  to  DCENTR 

The  inputs  for  base-jet  plume  interaction  effects  are  defined  using 
Option  1.  Incremental  forces  and  moments  due  to  jet  induced  boattail 
separation  and  separation  locations  on  aft  fins  are  calculated  if  these  inputs 
are  used. 


•  This  option  should  only  be  run  for  supersonic  cases 
(i.e.  Moo  ^  1.2) 

•  The  calculations  will  be  done  for  three  types  of  aft  bodies 
conical  boattail,  ogival  boattail,  or  cylindrical  (i.e.  no  boattail). 
Error  messages  will  be  printed  to  the  output  file  and  the 
calculations  skipped  if  any  other  aft  body  is  defined. 

•  If  BASE=. FALSE,  or  is  not  input  the  calculations  will  be 
skipped. 

•  DEXIT  must  not  equal  zero  if  this  option  is  used. 
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•  The  jet  Mach  number  (JMACH),  jet  to  freestream  static 
pressure  ratio  (PRAT),  and  jet  to  freestream  stagnation 
temperature  ratio  (TRAT)  must  be  specified  for  each 
freestream  Mach  number  or  velocity  input  in  the  namelist 
FLTCON.  Fo.  subsonic  or  transonic  freestream  Mach 
numbers  or  velocities,  dummy  values  must  be  input  for 
JMACH,  PRAT,  and  TRAT.  The  user  must  be  careful  to 
match  these  inputs  with  the  proper  freestream  conditions. 

•  If  a  portion  of  the  fins  in  a  fin  set  are  located  on  the  boattail  or 
base,  the  boattail  separation  locations  will  be  calculated  and 
output  at  each  fin  roll  angle.  However,  if  the  fins  do  not 
extend  to  the  boattail  the  separation  locations  will  be  skipped. 

•  Results  may  be  inaccurate  if  excessive  extrapolation  is 
required.  If  extrapolation  occurs,  a  warning  message  will  be 
printed  to  the  output  file.  To  avoid  extrapolation  and 
minimize  inaccuracy,  the  input  parameters  should  be  kept 
within  the  ranges  shown  in  Figure  5. 

314  NAMELIST  ELLBOD  -  Elliptical  Body  Geometry 

Elliptically-shaped  cross  section  bodies  are  defined  using  this  namelist. 
The  inputs  are  similar  to  those  for  the  axisymmetric  body  geometry 
(AXIBOD),  and  are  shown  in  Figures  6a  and  6b.  The  types  of  shapes  available, 
and  the  limitations,  are  the  same  as  those  given  for  axisymmetric  bodies. 
However,  the  base-jet  plume  interaction  input  options  in  namelist  AXIBOD 
are  not  available  in  namelist  ELLBOD.  Please  read  Section  3.1.3  for 
limitations. 

Note  that  the  body  cross  section  ellipticity  can  vary  along  the  body 
longitudinal  axis  Sections  which  are  taller-than-wide  and  wider-than-tall  can 
be  mixed  to  produce  "shaped"  designs.  The  shap)e  of  the  sections  is  controlled 
by  the  variables  ENOSE,  ECENTR,  and  EAFT  or  ELLIP,  H  and  W. 

3.1.5  NAMELIST  PROTUB  -  Protuberance  Geometry 

Missile  protuberances  can  be  input  using  this  namelist.  Axial  force 
coefficient  is  calculated  for  the  protuberances  and  added  to  the  body  axial  force 
coefficient.  Figure  7  shows  the  inputs  required.  Figure  8  shows  the  different 
protuberance  shapes  available.  Tlie  following  defines  the  inputs  required  for 
protuberance  calculations: 

•  NPROT  is  the  number  of  protuberance  sets.  A  protuberance 
set  is  made  up  of  protuberances  at  the  same  axial  location 
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with  the  same  size  and  shape.  Therefore,  it  is  only  necessary 
to  describe  the  geometry  of  one  individual  protuberance  per 
set.  The  maximum  number  of  protuberance  sets  is  20. 

•  NLOC  is  the  number  of  protuberances  in  each  protuberance 
set.  NLOC  accounts  for  the  number  of  identical  protuberances 
located  around  the  missile  body  at  a  given  axial  location. 

•  The  following  equation  helps  to  clarify  the  relationship 
between  NLOC  and  NPROT: 

NLOC(l)+NLOC(2)-;-NLOC(3)+  . . .  +NLOC(NPROT)  =  (Total 
number  of  protuberances  on  the  missile) 

•  The  axial  location  of  a  protuberance  (XPROT)  should  be  input 
at  the  protuberance  geometric  centroid.  An  approximation  of 
the  centroid  will  be  adequate  for  the  analysis.  The  location  is 
used  to  calculate  the  average  boundary  layer  thickness  over 
the  protuberance  length. 

•  VCYL,  HCYL,  BLOCK,  and  FAIRING  type  protuberances  have 
1  member.  LUG  types  have  4  members  and  SHOE  types  have 
3  members.  (Refer  to  Figure  8) 

•  All  inputs  for  LPROT,  WPROT,  HPROT,  and  OPROT  are  in 
sequential  order  based  upon  the  members  specified  with  the 
protuberance  type  (PTYPE)  input. 

•  The  FAIRING  type  protuberance  should  always  have  a  zero 
offset.  The  code  will  assume  a  zero  offset  even  if  a  non-zero 
offset  is  input. 

More  complex  protuberance  shapes  can  be  analyzed  by  a  component 
build-up  method.  Each  member  is  treated  as  a  separate  protuberance. 
Combinations  of  vertical  cylinders,  horizontal  cylinders,  and  flat  plates  or 
blocks  can  be  input  at  specified  offsets  from  the  missile  body.  If  a  FAIRING 
type  protuberance  is  used  in  a  component  build-up,  the  offset  should  be  zero. 
The  user  must  manually  add  axial  force  of  the  individual  members  of  the 
component  build-up  if  the  total  protuberance  axial  force  is  desired. 

Figure  9  shows  an  example  input  file  for  a  missile  with  several 
protuberance. 

3.1.6  NAMELIST  FINSETn  -  Define  Fin  Set  n 
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Figure  10a  describes  the  variables  needed  to  be  input  for  fin  set 
planform  geometry  descriptions.  Optional  fin  cross-section  inputs  are 
described  in  Figure  10b.  Special  user  specified  fin  cross-sections  can  be  input 
using  the  variables  in  Figure  10c.  The  user  may  specify  up  to  four  non¬ 
overlapping  fin  sets.  The  variable  "n"  in  the  namelist  specifies  the  fin  set 
number.  Fin  sets  must  be  numbered  sequentially  from  the  front  to  the  back  of 
the  missile  beginning  with  fin  set  one.  An  input  error  will  occur  if  "n"  is  zero 
or  omitted.  The  code  allows  for  between  1  and  8  geometrically  identical  panels 
to  be  input  per  fin  set.  The  panels  may  be  arbitrarily  rolled  about  the  body  and 
can  be  given  dihedral. 

Four  types  of  airfoil  sections  are  permitted—hexagonal  (HEX),  circular 
arc  (ARC),  NACA  airfoils  (NACA),  and  user  defined  (USER).  Only  one  type  of 
airfoil  section  can  be  specified  per  fin  set,  and  this  type  is  used  for  all  chord 
wise  cross  sections  from  root  to  tip.  Diamond-shaped  sections  are  considered  a 
special  case  of  the  HEX  type;  hence,  hexagonal  and  diamond  sections  can 
coexist  on  the  same  panel.  The  airfoil  proportions  can  be  varied  from  span 
station  to  span  station. 

The  user  selects  "break  points’"  on  the  panel  (Figure  11).  A  "break 
point"  specifies  a  change  in  leading  or  trailing  edge  sweep  angle.  Also  a  break 
point  may  specify  a  change  in  airfoil  section,  but  the  section  must  be  of  the 
same  type  (i.e.,  a  change  in  section  type  cannot  go  from  a  NACA  to  an  ARC) 
only  the  proportions  can  change.  The  location  of  each  "break  point"  is  defined 
by  specifying  its  semi-span  station  (SSPAN)  from  the  vehicle  centerline  and 
distance  from  the  first  body  station  to  the  chord  leading  edge  (XLE).  The 
"break  point"  chord  leading  edge  array  (XLE)  can  be  defined  by  simply 
specifying  the  root  chord  leading  edge  [XLE(l)]  and  the  sweep  angles  of  each 
successive  panel  segment  if  the  semi-span  stations  are  input.  Note  that  only 
those  variables  that  uniquely  define  the  fin  need  to  be  entered.  Redundant 
inputs  can  lead  to  numerical  inconsistencies  and  subsequent  computational 
errors. 


The  panel  sweep  angle  (SWEEP)  can  be  specified  at  any  span  station  for 
each  segment  of  the  panels.  If  STA=0.,  the  sweep  angle  input  is  measured  at 
the  segment  leading  edge;  if  STA=1.,  the  sweep  angle  input  is  measured  at  the 
segment  trailing  edge.  Note  that  some  aerodynamic  methods  are  very 
sensitive  to  panel  sweep  angle.  For  small  span  fins,  small  errors  in  the 
planform  inputs  can  create  large  sweep  angle  calculation  errors.  It  is 
recommended  that  exact  sweep  angles  be  specified  wherever  possible;  for 
example,  if  the  panel  trailing  edge  is  unswept,  specifying  SWEEP=0.  and 
STA=1.  will  minimize  calculation  error.  Then  the  leading  edge  sweep  will  be 
computed  by  the  code  internally  using  the  SSPAN  and  CHORD  inputs. 

Since  all  panels  are  assumed  to  be  planar  (i.e.,  no  tip  dihedral),  all 
inputs  must  be  "true  view".  Once  the  planform  of  a  single  panel  is  defined. 
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all  fins  of  the  set  are  assumed  to  be  identical.  The  number  of  panels  present  is 
defined  using  the  variable  NPANEL.  Each  panel  may  be  rolled  to  an  arbitrary 
position  around  the  body  using  the  variable  PHIF.  PHIF  is  measured 
clockwise  from  top  vertical  center  (looking  forward  from  behind  the  missile) 
as  shown  in  Figure  12.  Each  panel  may  also  contain  a  constant  dihedral.  A 
panel  has  zero  dihedral  when  it  is  aligned  along  a  radial  ray  from  the 
centerline  (see  Figure  12).  The  variable  used  to  specify  dihedral  is  GAM.  GAM 
is  positive  if  the  panel  tip  chord  is  rotated  clockwise. 

Different  aerodynamics  will  be  computed  depending  upon  whether  the 
FLTCON  namelist  variable  PHI,  or  the  FINSETn  namelist  variable  PHIF,  is 
used  to  roll  the  geometry.  Figure  13  depicts  the  usage  of  the  roll  options.  The 
variable  "PHI"  means  that  the  body  axes  system  is  to  be  rolled  with  the 
missile  body,  whereas  PHIF  keeps  the  aerodynamics  in  a  non-rolled  body  axis, 
but  rather  locates  the  fin  positions  around  the  body.  PHIF  must  be  input  for 
each  panel,  while  PHI  rolls  the  whole  configuration. 

It  is  the  user’s  responsibility  to  assure  that  the  fins  are  (1)  on  the  body 
surface,  and  (2)  do  not  lie  internal  to  the  body  mold  line.  The  program  does 
not  check  for  these  peculiarities.  If  SSPAN(1)=0  is  input,  the  program  will 
assume  that  the  panel  semi-span  data  relative  to  its  root  chord  are  supplied. 
The  code  will  automatically  interpolate  the  body  geometry  to  place  the  panel 
on  the  body  surface  with  the  root  chord  parallel  to  the  body  centerline.  See 
Section  3.4  for  modeling  fins  on  body  segments  of  varying  radii. 

When  defining  more  than  one  fin  set,  the  fin  sets  must  never  have 
their  planforms  overlap  one  another.  There  must  be  sufficient  space  between 
the  forward  fin  trailing  edge  and  aft  fin  leading  edge  to  avoid  violating  the 
assumptions  made  by  the  aerodynamic  computations.  It  is  assumed  by  the 
aerodynamic  model  that  the  vortices  are  fully  rolled  up  w’hen  they  pass  the 
control  points  of  the  next  downstream  set  of  fins.  In  reality  the  vortex  sheet 
does  not  fully  roll  up  until  it  is  at  least  four  semispans  downstream.  If  two  fin 
sets  are  closer  than  this  the  results  may  be  in  error  since  the  use  of  a  vortex 
filament  model  may  introduce  too  much  vorticity.  The  closer  the  spacing  the 
larger  the  error  may  be.  No  algorithm  error  will  result  from  too  close  a  fin  set 
spacing. 

Panels  with  cut-out  portions  can  be  modeled  by  using  one  of  the  ten 
available  fin  segments  as  a  transition  segment  This  is  accomplished  by  giving 
the  segment  a  small  span,  such  as  0.(X)C1,  and  specifying  the  segment  root  and 
tip  chords  to  transition  into  the  cut-out  portion  of  the  fin. 

3.1.7  NAMELIST  DEFECT  -  Panel  Deflection  Angles 

This  namelist  permits  the  user  to  fix  the  incidence  angle  for  each  panel 
in  each  fin  set.  The  variables  are  given  in  Figure  14.  Note  that  the  panel 
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numbering  scheme  is  assumed  to  be  that  shown  in  Figure  12.  The  array 
element  of  each  deflection  array  corresponds  to  the  panel  number. 

The  scheme  for  specifying  deflection  angles  is  uirique,  yet  concise.  The 
scheme  used  is  based  upon  the  body  axis  rolling  moment: 

"In  Missile  Datcom  a  positive  panel  deflection  is  one  which  will 
produce  a  negative  (counterclockwise  when  viewed  from  the 
rear)  roll  moment  increment  at  zero  angle  of  attack  and 
sideslip." 

3.1.8  NAMELIST  TRIM  -  Trim  Aerodynamics 

This  namelist  instructs  the  program  to  statically  trim  the  vehicle 
longitudinally  (Cm=0)  The  inputs  are  given  in  Figure  15.  Note  that  only  one 
fin  set  can  be  used  for  trimming.  The  user  only  specifies  the  range  of 
deflection  angles  desireii  using  DELMIN  and  DELMAX;  the  code  will  try  to 
trim  the  vehicle  for  each  angle  of  attack  specified  using  the  allowable  fin 
deflections.  This  option  will  not  trim  the  vehicle  a  specific  angle  of  attack  if 
the  deflection  required  is  outside  the  range  set  by  the  values  of  DELMIN  and 
DELMAX. 

The  deflection  sign  convention  used  is  that  described  in  Section  3.1.7; 
hence,  DELMIN  and  DELMAX  are  input  as  if  deflecting  the  panel  to  the 
maximum  will  produce  a  negative  rolling  moment  from  the  panels  resulting 
normal  force  increment.  DELMIN  must  always  be  less  than  DELMAX. 

A  logical  variable,  ASYM,  has  been  included  to  permit  reverse  panel 
deflections.  For  example,  deflecting  all  panels  in  one  sense  results  in  a  rolling 
moment  and  no  pitching  moment.  The  ASYM  flag  will  permit  analysis  of  an 
elevator  (or  pitch  deflection)  effect,  by  deflecting  panels  on  one  side  of  the 
vehicle  only,  with  opposite  panels  mirroring  those  deflections.  Since  a 
maximum  of  eight  panels  are  allowed  in  a  fin  set,  only  four  panels  of  the  fin 
set  can  be  deflected  in  the  reverse  direction  using  the  ASYM  flag.  Both 
trimmed  and  untrimmed  results  are  available  for  output. 

3.1.9  NAMELIST  INLET  -  Axisymmetric  and  2- Dimensional  Inlet  Geometry 

This  namelist  is  used  to  model  the  inlet  and  diverter  geometry. 
Axisymmetric,  two-dimensional  side  mounted,  and  two-dimensional  top 
mounted  external  compression  inlets  can  be  described.  The  inlets  may  be 
covered  or  uncovered  and  oriented  in  any  position  about  the  missile  body. 
Inlet  normal  force,  pitching  moment,  side  force,  yawing  moment,  and  axial 
force  are  calculated.  The  methods  are  valid  for  subsonic,  transonic,  and 
supersonic  speeds  Figure  16  shows  the  INLET  namelist  inputs,  and  Figures 
17a,  17b,  and  17c  show  the  inlet/ diverter  geometry  for  each  type  of  inlet. 
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•  Inlet  roll  orientation  uses  the  same  convention  as  the  fin 
panel  roll  orientation. 

•  Inlet  height  and  width  or  inlet  diameter  is  input  at  five  axial 
locations  desaibed  in  Figures  17a,  17b,  and  17c: 

1)  leading  edge  or  tip 

2)  cowl  lip  leading  edge 

3)  midbody  start 

4)  boattail  start 

5)  boattail  end 

•  If  the  inlet  is  covered  (COVER=.TRUE.),  no  flow  is  allowed 
into  the  inlet.  The  inlet  is  plugged  between  stations  1  and  2, 
flush  with  the  inlet  face. 

Inlet  additive  drag  or  spillage  drag  can  be  calculated  for  external 
compression  inlets  operating  at  off-design  conditions  (Moo<Mdesign)  for  Mach 
numbers  greater  than  1.  Whenever  flow  spillage  occurs,  the  mass  flow  ratio  is 
less  than  one,  and  additive  forces  are  generated  on  the  deflected  streamtube 
captured  by  the  inlet.  If  the  inlet  operates  on-design,  the  ramp  shock  lies  on 
the  inlet  face  and  on  the  cowl  lip.  In  these  cases,  the  maximum  mass  flow 
ratio  is  one  (zero  spillage)  and  the  minimum  additive  forces  are  zero. 

•  If  the  inlet  is  covered  (COVER=.TRUE.),  the  additive  drag 
calculations  will  be  skipped. 

•  If  ADD=. FALSE.,  or  is  not  input  the  additive  drag  calculations 
will  be  skipped. 

•  Mass  flow  ratio  (MFR)  must  be  specified  for  each  freestream 
Mach  number  or  velocity  given  in  namelist  FLTCON.  For 
Mach  numbers  less  than  1,  dummy  values  must  be  input  for 
MFR.  The  user  must  be  careful  to  match  these  inputs  with 
tl^e  proper  freestream  conditions. 

•  The  additive  drag  is  calculated  at  zero  angle  of  attack  and 
assumed  to  remain  constant  for  all  angles  of  attack. 

3.1.10  NAMELIST  EXPR  -  Experimental  Data  Substitution 

This  namelist  is  used  to  substitute  experimental  data  for  the  theoretical 
data  generated  by  the  program.  The  variables  to  be  input  are  shown  in  Figure 
18.  Use  of  namelist  EXPR  does  not  stop  the  program  from  calculating 
theoretical  data,  but  rather  the  experimental  data  is  used  in  configuration 
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synthesis,  and  it  is  the  experimental  data  that  is  used  for  the  component 
aerodymamics  for  which  it  is  input. 

Experimental  data  may  be  substituted  for  any  configuration  component 
or  partial  configuration.  Experimental  data  is  input  at  a  specific  Mach 
number,  when  using  namelist  EXPR  the  case  must  be  run  at  the  Mach 
number  for  which  you  are  substituting  experimental  data.  However,  the 
experimental  data  being  input  may  have  different  reference  quantities  and  a 
different  center  of  gravity  location  than  the  case  being  nm. 

Experimental  data  input  for  a  fin  alone  is  input  as  panel  data,  not  as 
total  fin  set  data.  The  user  should  note  that  experimental  data  for  fin  alone 
Cmo  is  not  used  in  the  configuration  synthesis  process.  Instead  fin  alone  Cno 
(the  experimental  value  if  input)  is  used  to  determine  the  fin  contribution  to 
Cma  during  configuration  synthesis.  If  body  alone  experimental  data  and  body- 
fin  experimental  data  are  input  for  the  same  case  the  body  data  is  ignored  in 
configuration  synthesis.  If  experimental  Cma  data  is  input  for  a  body  +  1  fin  set 
for  a  multi-fin  set  configuration,  the  calculated  contributions  to  Cma  of  the 
other  fin  sets  are  added  to  the  experimental  data. 

Since  the  experimental  namelist  forms  the  basis  for  configuration 
incrementing,  the  lateral  directional  coefficients  are  included  to  allow  for 
sideslip  cases.  These  coefficients  are  input  the  same  as  the  longitudinal 
coefficients.  However,  if  the  lateral  directional  coefficients  are  input,  the 
lateral  directional  beta  derivatives  will  not  be  computed  our  output. 

The  following  rules  apply  to  the  use  of  namelist  EXPR. 

•  It  is  assumed  that  the  coefficients  in  EXPR  are  for  the  same 
sideslip  and/or  aerodynamic  roll  as  the  case  being  run. 

•  Separate  namelist  EXPR  must  be  specified  for  each  Mach 
number. 

•  Each  namelist  EXPR  must  end  with  a  $END  card. 

•  Separate  namelist  EXPR  must  be  specified  for  each  partial 
configuration  for  which  experimental  data  is  to  be  input,  (i.e., 
body,  body  + 1  fin  set,  etc) 

•  Separate  namelist  EXPR  must  be  specified  for  each  reference 
quantity  change. 

Example: 

The  user  has  experimental  data  available  for  a  body  +  2  fin  set 
configurations  and  is  interested  in  the  effects  of  adding  a 
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booster  containing  a  third  fin  set.  he  would  then  use  namelist 
EXPR  to  input  the  experimental  data.  When  the 
configuration  is  synthesized,  it  would  use  the  experimental 
data  for  body  +  2  fin  sets  and  theoretical  data  for  fin  set  three. 

3.2  Control  Card  Inputs 

Control  cards  are  one  line  con  mands  which  select  program  options. 
Although  they  are  not  required  inputs,  they  permit  user  control  over 
program  execution  and  the  types  of  output  desired.  Control  cards  enable  the 
following: 

•  Printing  internal  data  array  results  for  diagnostic  purposes 
(DUMP) 

•  Outputting  intermediate  calculations  (PART,  BUILD, 
PRESSURES,  PRINT  AERO,  PRINT  EXTRAP,  PRINT  GEOM, 

PLOT,  NAMELIST,  WRITE,  FORMAT) 

•  Selecting  the  system  of  units  to  be  used  (DIM,  DERIV) 

•  Defining  multiple  cases,  permitting  the  reuse  of  previously 
input  namelist  data  or  deleting  namelists  of  a  prior  case 
(SAVE,  DELETE,  NEXT  CASE) 

•  Adding  case  titles  or  comments  to  the  input  file  and  output 
pages  (^  CASEID) 

•  Limits  the  calculations  to  longitudinal  aerodynamics  (NO 
LAT) 

3.2.1  Control  Card  -  General  Remarks 


A  total  of  42  different  control  cards  are  available.  There  is  no  limit  to 
the  number  of  control  cards  that  can  be  present  in  a  case.  If  two  or  more 
control  cards  contradict  each  other,  the  last  control  card  input  will  take 
precedence.  All  control  cards  must  be  input  as  shown,  including  any  blanks. 
Control  cards  can  start  in  any  column  but  they  cannot  be  continued  to  a 
second  card.  Misspelled  cards  are  ignored.  Control  cards  can  be  located 
anywhere  v.’ithin  a  case. 

Once  input,  the  following  control  cards  remain  in  effect  for  all 
subsequent  cases: 

DIM  FT  DIM  IN  DIM  CM  DIMM  FORMAT 

HYPER  INCREMENT  NOGO  NO  LAT  PLOT 
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SOSE 


WRITE 


The  following  control  cards  are  effective  only  for  the  case  in  which  they 
appear; 


BUILD 
DUMP  CASE 
PRESSURES 
PRINT  GEOM 


CASEID 
DUMP  NAME 
PRINT  AERO 
SAVE 


DAMP  DELETE 

NAMELIST  PART 

PRINT  EXTRAP 
SPIN  TRIM 


These  control  cards  can  be  changed  fron\  case  to  case: 


DERTV  DEG  DERIV  RAD  NACA 


The  only  control  card  that  can  be  optionally  saved,  from  case-to-case,  is 
the  NACA  card. 


3.2.2  Control  Card  Definition 


Available  control  cards  are  sunamiirized  as  follows: 

BUILD 

This  control  card  instructs  the  program  to  print  the  results  of  a 
configuration  build-up.  All  configurations  which  can  be  built  from  the 
components  defined  will  be  synthesized  and  output,  including  isolated  data 
(e.g.,  body  alone,  fin  alone,  etc.).  Component  build-up  data  is  not  provided  if 
the  TRIM  option  is  selected. 

CASEID 


A  user  supplied  title  to  be  printed  on  each  output  page  is  specified.  Up 
to  72  characters  Ccin  be  specified  (card  columns  8  to  80). 

DAMP 


When  DAMP  control  card  is  input  longitudinal  dynamic  derivatives 
are  computed  and  the  results  output  for  the  configuration.  Dynamic 
derivatives  for  configuration  components  or  partial  configurations  may  be 
output  using  the  PART  or  BUILD  control  cards  respectively. 

DELETE  namel.nar.te2 

This  control  card  instructs  the  program  to  ignore  a  previous  case 
namelisi  input  that  was  retained  using  the  SAVE  control  card.  All  previously 
saved  namelists  with  the  names  specified  will  be  purged  from  the  input  file. 


23 


Any  new  inputs  of  the  same  namelist  will  be  retained.  At  least  one  name 
(namel)  must  be  specified. 

DIM  IN.  DIM  FT.  DIM  CM,  or  DIM  M 

This  control  card  sets  the  system  of  units  for  the  user  inputs  and 
program  outputs.  The  four  options  are  inches  (DIM  IN),  feet  (DIM  FT), 
centimeters  (DIM  CM),  and  meters  (DIM  M).  The  default  system  of  units  is 
feet.  Once  the  system  of  units  has  been  set,  it  remains  set  for  all  subsequent 
cases  of  the  "run". 

DERIV  DEG  or  DERIV  RAD 

All  output  derivatives  are  set  to  either  degree  (DERIV  DEG)  or  radian 
(DERIV  RAD)  measure.  The  default  setting  is  degree.  The  derivative  units 
can  be  changed  more  than  once  during  the  run  by  inputting  multiple  DERTV 
cards. 

DUMP  CASE 


Internal  data  blocks,  used  in  the  computation  of  the  case,  are  written  on 
Tape  6.  This  control  card  automatically  selects  partial  output  (PART). 

DUMP  namel.name2 

This  permits  the  user  to  write  selected  internal  data  blocks  or  common 
blocks  on  Tape  6.  At  least  one  name  (namel)  must  be  specified.  The  arrays 
will  be  dumped  in  units  of  feet,  pounds,  degrees  or  degrees  Rankine.  Table  5 
shows  the  common  block  dump  names. 

FORMAT  (format) 

This  control  card  is  used  in  conjunction  with  the  WRITE  control  card. 
It  specifies  the  format  of  the  data  to  be  printed  to  tape  unit  3.  The  format  is 
input  starting  with  a  left  parenthesis,  the  format  and  a  right  parenthesis.  This 
is  exactly  the  same  as  a  FORTRAN  FORMAT  statement.  Because  of  the  code 
structure,  alphanumeric  data  must  not  be  printed.  For  example; 

FORMAT  ((8(2X,F10.4))  is  legal 

FORMAT  (2HX=,F10.4)  is  illegal 

The  default  format  is  8F10.4,  and  will  be  used  if  the  FORMAT  control 
card  is  not  present.  Multiple  formats  can  be  used.  The  last  FORMAT  read  will 
be  used  for  all  successive  WRITE  statements  until  another  FORMAT  is 
encountered.  Hence,  the  FORMAT  must  precede  the  applicable  WRITE. 

HYPER 
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This  control  card  causes  the  program  to  select  the  Newtonian  flow 
method  for  bodies  at  any  Mach  number  above  1 .4.  HYPER  should  normally  be 
selected  at  Mach  niunbers  greater  than  6. 

INCRMT 

This  card  is  used  to  set  the  configuration  incrementing  flag. 
Configuration  incrementing  uses  the  first  case  of  a  run  to  determine 
correction  factors  for  the  longitudinal  and  lateral  aerodynamic  coefficients. 
These  correction  factors  are  computed  by  comparing  theoretical  and 
experimental  values  for  each  coefficient  for  which  data  is  input.  The 
experimental  values  are  input  using  namelist  EXPR.  During  subsequent  cases 
of  the  run,  the  correction  factors  are  applied  to  coefficients  for  which 
experimental  data  was  input  in  the  first  case.  This  provides  the  user  with  a 
method  to  evaluate  changes  in  a  configuration. 

The  INCRMT  card  must  be  input  in  the  first  case  of  a  run.  The  first  case 
must  be  run  at  the  same  Mach  number  as  the  experimental  data  which  is 
input.  Once  the  increment  flag  is  set  it  cannot  be  deleted  during  that  run. 

The  following  restrictions  apply: 

•  All  cases  of  a  run  must  have  the  same  number  of  fin  sets. 

•  All  cases  of  a  run  must  have  the  same  sideslip  or 
aerodynamic  roll  angle  as  the  first  case  (BETA  or  PHI  as 
specified  in  namelist  FLTCON). 

•  The  first  case  must  be  run  at  exactly  the  same  angles  of  attack 
as  the  experimental  data  being  input. 

•  All  cases  must  be  run  within  the  same  Mach  regime 
(subsonic,  transonic,  or  supersonic)  as  the  experimental  data. 

•  Experimental  data  can  only  be  input  in  the  first  case  and  only 
for  the  complete  configuration.  No  additional  data  can  be 
substituted. 

•  To  increment  Cy^  and  Cnq  experimental  data  must  be  input 
for  Cy  and  Cn- 

Use  of  configuration  incrementing  may  or  may  not  increase  the 
accuracy  of  the  results.  The  following  guidelines  will  produce  better  results 
when  using  configuration  incrementing: 
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•  The  user  may  run  difterent  angles  of  attack  in  each  case. 
However,  no  angle  of  attack  should  exceed  the  upper  or 
lower  limit  of  the  angles  of  attack  for  which  experimental 
data  was  input  in  the  first  case. 

•  Experimental  data  should  be  input  at  as  many  angles  of  attack 
as  possible. 

•  The  user  should  remember  that  the  effect  of  a  change  in 
Mach  number  from  case  to  case  is  not  corrected  by  inputting 
experimental  data  at  one  Mach  number  as  is  required. 

NACA 

This  card  defines  the  NACA  airfoil  section  designation  (or  supersonic 
airfoil  definition).  Note  that  if  airfoil  coordinates  and  the  NACA  card  are 
specified  for  the  same  aerodynamic  surface,  the  airfoil  coordinate  specification 
will  be  used.  Therefore,  if  coordinates  have  been  specified  in  a  previous  case 
and  the  SAVE  option  is  in  effect,  the  saved  namelist  must  be  deleted  or  the 
namelist  variable  SECTYP  must  be  changed  for  the  NACA  card  to  be 
recognized  for  that  aerodynamic  surface.  The  airfoil  designated  with  this  card 
will  be  used  for  all  segments  and  panels  of  the  fin  set. 

The  form  of  this  control  card  and  the  required  parameters  are  as 
follows: 


Card  Column(s)  Input(s) 

Purpose 

1  thru  4 

NACA 

The  unique  letters  NACA 
designate  that  an  airfoil  is  to  be 
defined 

5 

Any  delimiter 

6 

1,2,3,  or  4 

Fin  set  number  for  which  the 
airfoil  designation  applies 

7 

Any  delimiter 

8 

1,4,5,6,S 

Type  of  NACA  airfoil  section; 
1-series  (1),  4-digit  (4),  5-digit  (5), 
6-series  (6),  or  supersonic  (S) 

9 

Any  delimiter 

10  thru  80 

Designation 

Input  designation  (see  Table  6); 
columns  are  free-field  (blanks 
are  ignored) 

Only  fifteen  (15)  characters  are  accepted  in  the  airfoil  designation.  The 
vocabulary  consists  of  the  following  characters: 

01234  5  6789A,= 
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Any  characters  input  that  are  not  in  the  vocabulary  list  will  be 
interpreted  as  the  number  zero  (0).  Table  6  details  the  restrictions  on  the 
NACA  designation. 

NAMELIST 

This  control  card  instructs  the  program  to  print  all  namelist  data.  This 
is  useful  when  multiple  inputs  of  the  same  variable  or  namelist  are  used. 

NEXT  CASE 


This  card  indicates  termination  of  the  case  input  data  and  instructs  the 
program  to  begin  case  execution.  It  is  required  for  multiple  case  "rims".  This 
card  must  be  the  last  card  input  for  the  case. 

NQGO 


This  control  card  permits  the  program  to  cycle  through  all  of  the  input 
cases  without  computing  configuration  aerodynamics.  It  can  be  present 
anywhere  in  the  input  stream  and  only  needs  to  appear  once.  This  option  is 
useful  for  performing  error  checking  to  insure  all  cases  have  been  correctly  set 
up. 

NO  LAT 

This  control  card  inhibits  the  calculation  of  the  lateral-directional 
derivatives  due  to  sideslip  angle.  Savings  in  computation  time  can  be  realized 
by  using  this  option.  This  option  is  automatically  selected  when  using  TRIM. 

PART 


This  control  card  permits  printing  of  partial  aerodynamic  output,  such 
as  a  summary  of  the  normal  force  and  axial  force  contributors.  Partial  output 
of  the  configuration  synthesis  methods  is  only  provided  if  the  TRIM  option  is 
not  selected.  Use  of  this  card  is  equivalent  to  inputing  all  PRINT  AERO  and 
PRINT  GEOM  control  cards  (except  +). 

PLOT 


A  data  file  for  use  with  a  post-processing  plotting  program  is  provided 
when  this  control  card  is  used.  A  formatted  file  is  written  to  unit  3.  Appendix 
B  shows  the  format  of  this  data  file. 
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PRESSURES 


This  control  card  instructs  the  program  to  print  the  body  and  fin  alone 
pressure  coefficient  distributions  at  supersonic  speeds.  Only  pressure  data  to 
15  degrees  angle  of  attack  for  bodies  and  at  zero  angle  of  attack  for  fins  are 
printed. 

PRINT  AERO  name 


This  control  card  instructs  the  program  to  print  the  incremental 
aerodynamics  for  "name ",  which  can  be  one  of  the  following; 


BODY 

for  body  aerodynamics 

FINl 

for  FINSETl  aerodynamics 

FIN2 

for  FINSET2  aerodynamics 

FIN  3 

for  FINSET3  aerodynamics 

FIN4 

for  FINSET4  aerodynamics 

SYNTHS 

for  configuration  synthesis  aerodynamics 

TRIM 

for  trim /untrimmed  aerodynamics 

BEND 

for  panel  bending  moments 

HINGE 

for  panel  hinge  moments 

INLET 

for  inlet  aerodynamics 

All  options  are  automatically  selected  when  the  control  card  PART  is 
used.  Details  of  the  output  obtained  v/ith  these  options  are  presented  in 
Section  4.2. 

PRINT  EXTRAP 


This  control  card  enables  the  printing  of  method  extrapolation 
messages  produced  during  execution  of  the  case.  Extrapolation  messages  are 
not  normally  provided. 

PRINT  GEOM  name 


This  control  card  instructs  the  program  to  print  the  geometric 
characteristics  of  the  configuration  component  "name",  which  can  be  one  of 
the  following: 

BODY  for  body  geometry 

FINl  for  FINSETl  geometry 

F1N2  for  FINSET2  geometry 

FIN3  for  FINSET3  geometry 

F1N4  for  FINSET4  geometry 

INLET  for  inlet  geometry 
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All  options  are  automatically  selected  when  the  control  card  PART  is 

used. 

SAVE 


The  SAVE  card  saves  namelist  inputs  from  one  case  to  the  following 
case  but  not  for  the  entire  run.  This  permits  the  user  to  build-up  or  change  a 
complex  configuration,  case-to-case,  by  adding  new  namelist  cards  without 
having  to  re-input  namelist  cards  of  the  previous  case.  When  changing  a 
namelist  that  has  been  saved,  the  namelist  must  first  be  deleted  using  the 
delete  control  card. 

The  only  control  card  that  can  be  optionally  saved,  case-to-case,  is  the 
NACA  card. 

SOSE 


The  presence  of  this  control  card  selects  the  Second-Order  Shock 
Expansion  Method  for  axisymmetric  bodies  at  supersonic  speeds.  SOSE 
should  be  selected  if  any  Mach  number  is  higher  than  2.0. 

SPIN 


When  the  SPIN  control  card  is  input,  spin  and  magnus  derivatives  are 
computed  for  body  alone.  If  the  configuration  being  run  is  a  body  +  fin  sets, 
the  spin  derivatives  are  still  computed  for  body  alone.  A  PART  or  BUILD  card 
must  be  input  for  body  alone  derivatives  to  be  printed  out. 

TRIM 


This  control  card  causes  the  program  to  perform  a  trim  calculation. 
Component  buildup  data  cannot  be  dumped  if  TTUM  is  selected.  The  use  of 
this  control  card  is  the  same  as  if  the  namelist  TRIM  was  included  except  that 
the  defaults  for  namelist  TRIM  are  used. 

WRITE  name,  start.end 

This  control  card  causes  the  common  block  ’hame"  to  be  printed  to 
tape  unit  3  using  the  most  recent  FORMAT  control  card.  Locations  from 
"start”  to  "end"  are  dumped  (see  Table  5  for  common  block  write  names). 
Multiple  WRITE  statements  may  be  input,  and  there  is  no  limit  to  the 
number  which  may  be  present  The  presence  of  a  WRITE  will  cause  the  block 
"name"  to  be  printed  for  all  cases  of  the  run.  The  output  will  be  in  units  of 
feet,  pounds,  degrees,  or  degrees  Rankine.  If  the  PLOT  option  is  also  selected, 
this  output  will  be  "mixed"  with  the  PLOT  file  output  on  tape  unit  3. 
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Any  card  with  an  asterisk  (*)  in  Column  1  will  be  interpreted  as  a 
comment  card.  This  permits  detailed  documentation  of  case  inputs. 

3.3  Typical  Case  Set-up 

Figure  19  schematically  shows  how  Missile  Datcom  inputs  are 
structured.  This  example  illustrates  a  multiple  case  job  in  which  case  2  uses 
part  of  the  case  1  inputs.  This  is  accomplished  through  use  of  the  SAVE 
control  card.  Case  1  is  a  body-wing-tail  configuration;  partial  output, 
component  buildup  data,  and  a  plot  file  are  created.  Case  2  uses  the  body  and 
tail  data  of  case  1  (the  wing  is  deleted  using  DELETE),  specifies  panel 
deflection  angles  and  sets  the  data  required  to  trim. 

There  is  no  limit  to  the  number  of  cases  that  can  be  "stacked "  in  a 
single  run,  provided  that  no  more  than  300  namelist  inputs  are  "saved" 
between  cases.  If  a  SAVE  control  card  is  not  present  in  a  case,  all  previous  case 
inputs  are  deleted. 

3.3. 1  Configuration  Incrementing  Case  Set-up 

A  "configuration  incrementing"  case  setup  is  shown  in  Figure  20.  This 
figure  shows  the  inputs  for  a  three  case  setup  fin  parametric  analysis.  The  first 
case  is  the  calibration  case  with  the  remaining  cases  being  used  for  the 
parametric  analysis.  ITierefore,  the  first  case  must  contain  both  the  INCRMT 
control  card  and  EXPR  namelist.  These  should  only  appear  in  the  first  case. 

3.4  Special  Usage  of  Input  Parameters 

It  is  possible  to  manipulate  the  input  geometry,  such  that  unique 
configurations  can  be  modeled.  This  section  defines  those  capabilities. 

3.4.1  Locating  Panels  on  Varying  Body  Radii  Segments 

The  fin  panels  can  be  located  anywhere  on  the  geometry.  If  they  are  to 
be  positioned  on  a  varying  radii  .segment,  select  the  root  chord  span  station 
(SSPAN(l)J  such  that  the  center  of  the  exposed  root  chord  is  on  the  surface 
mold  line.  Physically  this  places  part  of  the  panel  within  the  body  and  part 
offset  from  the  body. 

If  SSPAN(l)  is  input  precisely  as  zero,  the  code  will  assume  that  panel 
semi-span  stations  relative  to  the  root  chord  are  defined.  It  will  then 
interpolate  the  body  geometry  at  the  root  chord  center  and  add  the  body 
radius  at  this  point  to  the  user  defined  values  m  the  SSPAN  array. 


30 


Table  1  Body  Addressable  Configurations 


CONFIGURATION 

SUPERSONIC 

M>1.2 

1.  Nose  Shape 

Conical 

Sharp 

X 

X 

X 

Blunted 

X 

X 

X 

Truncated 

X 

X 

X 

Tangent  Ogive 

Sharp 

X 

X 

X 

Blunted 

X 

X 

X 

Truncated 

X 

X 

X 

Other 

X 

X 

2.  Centerbody  Shape 

Cylinder 

X 

X 

X 

Eliiptically  Variable 

X 

X 

X 

3.  Afterbody  Shape 

Boattails 

Conical 

X 

X 

X 

Tangent  Ogive 

X 

X 

X 

Other 

X 

Flares 

Conical 

X 

X 

X 

Ogive 

X 

Other 

X 
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Table  2  Subsonic/Transonic  Method  Limitations 


METHOD 

RANGE  PERMITTED 

ACTION  TAKEN 

IF  EXCEEDED 

Nose  Bluntness  (Cn,  Cm) 

Sharp  Only 

Uses  Sharp  Method 

Conical  Nose  Slope 

0  to  25  Degrees 

Uses  25  Degrees 

Boattail  Shape 

Cone  or  Ogive 

Uses  Cone 

Conical  Boattail  Slope 

0  to  16  Degrees 

Extrapolates 

Ogive  Boattail  Slope 

0  to  28  Degrees 

Extrapolates 

Flare  Shape 

Cone 

Uses  Cone 

Flare  Slope 

Oto  10  Degrees 

Extrapolates 

Airfoil  t/c 

0to12% 

Continues,  H  Possible 
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Table  3  Namelist  Alphanumeric  Constzmts 


NAMELIST 

PERMITTED 

ALPHANUMERIC 

CONSTANTS 

CONVERTED 

VALUE 

(ALL) 

UNUSED 

1.0E‘30  (Initialized  Value) 

REFQ 

TURB 

0. 

i 

NATURAL 

1. 

AXIBOD  or  ELLBOD 

CONICAL 

0. 

CONE 

0. 

OGIVE 

1. 

POWER 

2. 

HAACK 

3. 

KARMAN 

4. 

PROTUB 

VCYL 

1, 

HCYL 

2. 

BLOCK 

5. 

FAIRING 

6. 

LUG 

3. 

SHOE 

4. 

FINSETn 

HEX 

0. 

NACA 

1. 

ARC 

2. 

USER 

3. 

INLET 

AXI 

3. 

20TOP 

1. 

2DSIDE 

2. 

EXPR 

BODY 

1. 

1 

FI 

2. 

F2 

3. 

F3 

4. 

F4 

5. 

BF1 

6. 

BF12 

7. 

BF123 

8. 

BF1234 

9. 
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Table  4  Equivalent  Sand  Roughness 


TYPE  OF  SURFACE 

EQUIVALENT  SAND 
ROUGHNESS 
k  (INCHES) 

RHR 

Acrodynamically  Smooth 

0.0 

0.0 

Pollshod  fHIttal  or  Wood 

0.02  E‘3  to  O.OS  E*3 

6  to  26 

Natural  Shaet  Motal 

0.16  E-3 

53 

Smooth  Matta  Paint,  Carefully  Appllad 

0.25  E-3 

83 

StarKtard  Camouflage  Paint, 

Average  Application 

0.40  E-3 

133 

Camouflage  Paint, 
Maaa-Produetlon  Spray 

1.20  E-3 

400 

Dip  Galvanized  Metal  Surface 

6.0  E-3 

2000 

Natural  Surface  of  Caat  Iron 

10.0  E-3 

3333 

PREFERRED  RHR  VALUES 


APPLICATION 

RHR 

Steel  Structural  Parte 

250 

Aluminum  and  THanlum  Structural  Parte 

125 

Cloao  Tolerance  Surfaces 

63 

Seale 
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NAMELIST  FLTCON 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

DEFINITION 

UNITS 

DEFAULT 

NALPHA 

- 

NUMBER  OF  ANGLES  OF  ATTACK 
(AT  LEAST  2) 

- 

• 

ALPHA 

20 

ANGLE  OF  ATTACK  OR  TOTAL  ANGLE 

OF  ATTACK 

DEG 

* 

BETA 

- 

SIDESLIP  ANGLE 

DEG 

0. 

PHI 

- 

AERODYNAMIC  ROLL  ANGLE 

DEG 

0. 

NMACH 

- 

NUMBER  OF  MACH  NUMBERS  (AT  LEAST  1) 

- 

- 

MACH  I 

20 

MACH  NUMBERS 

- 

- 

REN 

20 

REYNOLDS  NUMBER  PER  UNIT  LENGTH 

■fife 

- 

ALT 

1,/T 

GEOMETRIC  ALTITUDE 

■80 

0. 

VINF 

In 

f  20 

FREESTREAM  SPEED 

L/SEC( 

9  ■ 

TINF 

20 

FREESTREAM  STATIC  TEMPERATURE 

DEG 

r  • 

PINF  J 

20 

FREESTREAM  STATIC  PRESSURE 

■■■■■A 

(P  THE  FOLLOWING  COMBINATIONS  SATISFY  THE  REYNOLDS  NUMBER  AND 
MACH  NUMBER  INPUT  REQUIREMENTS 


USER  INPUT  PROQRAM  COMPUTES 


1 .  MACH,  REN 

2.  MACH,  ALT 

3.  VINF,  ALT 

4.  VINF,  TINF,  PINF 

5.  MACH,  TINF,  PINF 


(NONE) 

PINF,  TINF,  REN 
PINF,  TINF.  MACH,  REN 
MACH.  REN 
VINF.  REN 


®  INPUT  AS  1/FT  FOR  ENGLISH  UNITS  AND  1/M  FOR  METRIC  UNITS 
(P  INPUT  AS  FT  FOR  ENGLISH  UNITS  AND  M  FOR  METRIC  UNITS 
0  INPUT  AS  FT/SEC  FOR  ENGLISH  UNITS  AND  M/SEC  FOR  METRIC  UNITS 
®  INPUT  AS  LB/FT^FOR  ENGLISH  UNITS  AND  N/M*  FOR  METRIC  UNITS 


Figure  1  Flight  Conditions  Inputs 
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NAMELIST  REFQ 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

DERNIHON 

UNITS 

DEFAULT 

SREF 

- 

REFERENCE  AREA 

L'L 

mpm 

LF  ;f 

- 

REFERENCE  LENGTH  (LONGITUDINAL) 

L 

LATREF 

- 

REFERENCE  LENGTH 

L 

LREF 

(LATERAL-DIRECTIONAL) 

roughL^ 

SURFACE  ROUGHNESS  HEIGHT 

L@ 

0. 

RHR  JS 

✓ 

ROUGHNESS  HEIGHT  RATING 

- 

0. 

XCG 

- 

LONGITUDINAL  POSITION  OF  C.G.  (+AFT) 

L 

0. 

ZCG 

- 

VERTICAL  POSITION  OF  C.G.  (+UP) 

L 

0. 

SCALE 

- 

VEHICLE  SCALE  FACTOR  (MULTIPLIER  TO 

1. 

GEOMETRY) 

BLAYER 

- 

BOUNDARY  LAYER  TYPE  TURB  FOR  FULLY 

. 

TURB 

TURBULENT.  NATURAL  FOR  NATURAL 

TRANSITION 

®  EITHER  CAN  BE  USED 

®  DEFAULT  IS  BODY  MAXIMUM  CROSS-SECTIONAL  AREA,  IF  NO  BODY  IS  INPUT, 
MAXIMUM  FIN  PANEL  AREA  IS  USED. 

®  DEFAULT  IS  BODY  MAXIMUM  DIAMETER.  IF  NO  BODY  IS  INPUT.  MAXIMUM  FIN  PANEL 
MEAN  GEOMETRIC  CHORD  IS  USED. 

@  INPUT  AS  INCHES  FOR  ENGLISH  UNITS  AND  CENTIMETERS  FOR  METRIC  UNITS, 


Figure  2  Reference  Quantity  Inputs 
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NAMELIST  AXIBOD 

OPTION  1  INPUTS 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

DERNIT10N 

UNITS 

DEFAULT 

XOORXO 

LONGITUDINAL  COORDINATE  AT 

L 

0. 

NOSE  TIP 

TNOSE 

NOSE  SHAPE  NAME;  CONICAL.  CONE, 

- 

OGIVE 

OGIVE,  POWER,  HAACK,  OR  KARMAN 

POWER 

EXPONENT,  n,  FOR  POWER  SERIES 

- 

0. 

SHAPES.(r/R)  =  (x/L)" 

LNOSE 

NOSE  LENGTH 

L 

- 

DNOSE 

NOSE  SECTION  BASE  DIAMETER 

L 

1.0 

BNOSE 

NOSE  BLUNTNESS  RADIUS  OR 

L 

0. 

TRUNCATED  NOSE  RADIUS 

TRUNC 

TRUNCATION  FLAG  (.TRUE.  IF  NOSE  IS 

- 

.FALSE. 

TRUNCATED) 

LCENTR 

CENTERBCDY  LENGTH 

L 

0. 

DCENTR 

CENTERBODY  BASE  DIAMETER 

L 

DNOSE 

TAFT 

AFTERBODY  SHAPE  NAME:  CONICAL. 

- 

CONICAL 

CONE.  OR  OGIVE 

LAFT 

AFTERBODY  LENGTH 

L 

0. 

DAFTj 

AFTERBODY  BASE  DIAMETER 

L 

- 

DEXIT 

DIAMETER  OF  NOZZLE  EXIT 

L 

* 

BASE 

BASE-JET  PLUME  INTERACTION  FLAG 

- 

.FALSE. 

(  TRUE.  IF  CALCULATIONS  DESIRED) 

BETAN 

i 

NOZZLE  EXIT  ANGLE 

DEG 

. 

jmach' 

20 

JET  MACH  NUMBER  AT  NOZZLE  EXIT 

- 

. 

PRAT 

m 

\  20 

JET  TO  FREESTREAM  STATIC  PRESSURE 

- 

■ 

RATIO 

TRAT 

4 

1 

20 

JET  TO  FREESTREAM  STAGNATION 

- 

- 

_ 

TEMPERATURE  RATIO 

CD  AFT  BODY  MUST  NOT  BE  CYLINDRICAL  (i.e.  DAFT  NOT  EQUAL  TO  DCENTR) 

(2)  ONLY  REQUIRED  IF  BASE-JET  PLUME  INTERACTION  CALCULATIONS  DESIRED, 
(DEXIT  MUST  NOT  EQUAL  ZERO) 

(3)  JMACH.  PRAT,  AND  TRAT  ARE  SPECIFIED  FOR  EACH  FREESTREAM  MACH 
NUMBER  OR  VELOCITY  INPUT  IN  NAMELIST  $FLTCON 


Figure  3a  Axisymmetric  Body  Geometry  Inputs  -  Option  1 
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NAMELIST  AXIBOD 

OPTION  2  INPUTS  (USE  ONLY  IF  MACH  >  1.2) 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

DEHNITK>N 

UNITS 

DEFAULT 

XO  OR  XO 

- 

LONGITUDINAL  COORDINATE  AT  NOSE  TIP 

0. 

BNOSE 

- 

NOSE  BLUNTNESS  RADIUS  OR  TRUNCATED 
NOSE  RADIUS 

H 

0. 

TRUNC 

- 

TRUNCATION  FLAG  (  TRUE.  IF  NOSE  IS 
TRUNCATED) 

■ 

FALSE. 

DEXIT 

- 

DIAMETER  OF  NOZZLE  EXIT 

- 

NX 

- 

NUMBER  OF  INPUT  STATIONS 

- 

X  O 

50 

LONGITUDINAL  COORDINATES 

- 

R 

50 

RADIUS  AT  EACH  X  STATION 

- 

DISCON 

20 

INDICES  OF  X  STATIONS  WHERE  THE 
SURFACE  SLOPE  IS  DISCONTINUOUS 

■ 

• 

(I)  X(NX)  MUST  BE  END  OF  BODY 


Figure  3b  Axisymmetric  Body  Geometry  Inputs  -  Option  2 
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—  -DNOSE- 


-DCENTO  ►+X 


XO-M4-*  LNOSE 


LCENTR • 


8NOSE 


BLUNTED  NOSE 


ONOSE 


BMOSE 


TRUNCATED  NOSE 


DNOSE 


LNOSE 


1^ — H 

' LNOSE ' 


DEXIT 


Figure  4  Body  Geometry  Inputs 
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Input  Pafamnter 

Symbol 

Min 

Value 

Max. 

Value 

Boattail  shape 

CylirKler.  Cone.  Ogive  I 

Boattall  fineness  falio 

170 

0 

2 

Boattail  terminal  angle 

0' 

12“ 

Jet  pressure  ratio 

Pj'P- 

0 

10 

Freostroam  Mach 
number 

2 

5 

Arrgifl  of  Attack 

a 

0’ 

8“ 

Jer  ivlacn  number 

Mj 

1 

M.^1 

Nozzle  lerminal  angle 

Pn 

5" 

25" 

Jet  diameter  ratio 

CU/Db 

0  80 

0.95 

Jet  temperature  ratio 

j  /T|(»o 

4 

10 

Note;  If  Input  parameter  Is  not  between  minimum  and 
maximum  value  the  code  will  extrapolate 


Figure  5  Base-Jet  Plume  Interaction  Parameters 
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NAMELIST  ELLBOD 

OPTION  1 1NPUTS 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

DEFINITION 

UNITS 

DEFAULT 

XO  OR  XO 

- 

l  ONGITUDINAL  COORDINATE  AT 

NOSE  TIP 

L 

0 

TNOSE 

• 

NOSE  SHAPE  NAME:  CONICAL,  CONE. 
OGIVE.  POWER,  HAACK.  OR  KARMAN 

- 

OGIVE 

POWER 

- 

EXPONENT,  n.  FOR  POWER  SERIES 
SHAPES.  (r/R)  -  (x^L)" 

* 

0. 

LNOSE 

- 

NOSE  LENGTH 

L 

WNOSE 

- 

NOSE  SECTION  BASE  WIDTH 

L 

1.0 

BNOSE 

NOSE  BLUNTNESS  RADIUS  OR 

TRUNCATED  NOSE  RADIUS 

L 

0. 

TRUNC 

* 

TRUNCATION  FLAG  (.TRUE.  IF  NOSE  IS 
TRUNCATED) 

• 

.FALSE. 

ENOSE 

- 

ELLIPTICITY  AT  NOSE  BASE  (H/W) 

- 

1.0 

LCENTR 

- 

CENTERBODY  LENGTH 

L 

0. 

WCENTR 

- 

CENTERBODY  BASE  WIDTH 

L 

WNOSE 

ECENTR 

- 

ELLIPTICITY  AT  CENTERBODY  BASE  (H/W) 

• 

1.0 

TAFT 

•G 

• 

AFTERBODY  SHAPE  NAME:  CONICAL. 
CONE,  OR  OGIVE 

• 

CONICAL 

LAFT 

- 

AFTERBODY  LENGTH 

L 

0. 

WAFT 

- 

AFTERBODY  BASE  WIDTH 

L 

. 

EAR  J 

- 

ELLIPTICITY  AT  AR  BODY  BASE  {HM) 

- 

1.0 

DEXIT 

- 

DIAMETER  OF  NOZZLE  EXIT 

L 

- 

a  AFT  BODY  MUST  NOT  BE  CYLINDRICAL  (i.a.  WAFT  NOT  EQUAL  TO  WCENTR) 


Figure  6a  Elliptical  Body  Geometry  Inputs  -  Option  1 
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NAMELIST  ELLBOD 

OPTION  2  INPUTS  (USE  ONLY  IF  MACH  >  1.2) 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

DEFINITION 

UNITS 

DEFAULT 

XOOR  XO 

LONGITUDINAL  COORDINATE  AT  NOSE  TIP 

L 

0. 

BNOSE 

- 

NOSE  BLUNTNESS  RADIUS  OR  TRUNCATED 
NOSE  RADIUS 

L 

0. 

TRUNC 

* 

TRUNCATION  FLAG  (  TRUE.  IF  NOSE  IS 
TRUNCATED) 

* 

FALSE. 

DEXIT 

- 

DIAMETER  OF  NOZZLE  EXIT 

L 

NX 

- 

NUMBER  OF  INPUT  STATIONS 

- 

X  Q 

50 

LONGITUDINAL  COORDINATES 

L 

W  (2) 

50 

BODY  HALF-WIDTH  AT  EACH  X  STATION 

L 

DISCON 

20 

INDICES  OF  X  STATIONS  WHERE  THE 
SURFACE  SLOPE  IS  DISCONTINUOUS 

“ 

H  (D^ 

50 

BODY  HALF-HEKSHT  AT  EACH  X  STATION 

L 

• 

ELLIP© 

50 

BODY  HEIGHT  TO  WIDTH  RATO  AT  EACH 

X  STATION 

• 

1.0 

0)  X{NX>  MUST  BE  END  OF  BODY 

®  ONE  OF  THE  FOLLOWING  COMBINATIONS  IS  REQUIRED: 
W  AND  H,  W  AND  ELLIP.  OR  H  AND  ELLIP 


Figure  6b  Elliptical  Body  Geometry  Inputs  •  Option  2 
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NAMELIST  PROTUB 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

DERNITKON 

UNITS 

DEFAUI'i 

NPR07 

- 

NUMBER  OF  PROTUBERANCE  GETS 
(20  MAXIMUM) 

- 

0. 

PTYPE 

20 

PROTUBERANCE  SET  TYPE:  VCYL,  HCYL, 
BLOCK,  PAIRING.  LUG.  OR  SHOE  © 

- 

■ 

XPROT 

20 

LONGITUDINAL  DISTANCE  FROM  MISSILE 
NOSE  TO  PROTUBERANCE  SET 

L 

• 

NLOC  © 

20 

NUMBER  OF  PROTUBERANCES  IN  EACH 
PROTUBERANCE  SET 

0. 

LPROT 

100 

LENGTH  OF  EACH  MEMBER  OR 
PROTUBERANCE 

L 

• 

WPROT 

-c 

O 

o 

/=\ 

WIDTH  OF  EACH  MEMBER  OR 
PROTUBERANCE 

L 

- 

HPROT 

100 

HEIGHT  OF  EACH  MEMBER  OR 
PROTUBERANCE 

L 

• 

OPROT  J 

100 

VERTICAL  OFFSET  OF  EACH  MEMBER 

OR  PROTUBERANCE 

L 

0. 

®  NLOC  ACCOUNTS  FOR  IDENTICAL  PROTUBERANCES  (SAME  SIZE  AND  SHAPE)  LOCATED 
AROUND  THE  BODY  AT  THE  SAME  AXIAL  LOCATION 
<2)  LUG  TYPE  HAS  4  MEMBERS.  SHOE  TYPE  HAS  3  MEMBERS.  LPROT,  WPROT,  HPROT.  AND 
OPROT  MUST  BE  SPECIFIED  FOR  EACH  MEMBER. 

0)  INPUT  FOR  EACH  PROTUBERANCE  (VCYL,  HCYL.  BLOCK.  OR  FAIRING  TYPE)  OR  EACH 
PROTUBERANCE  MEMBER  (LUG  AND  SHOE  TYPE) 


Figure  7  Protuberance  Inputs 


Vertical  Circular  Cylinder 
Oriented  Peipendicular  To 
Flow  (VCYL): 


Horizontal  Circular 
Cylinder  Oriented  In  Flow 
Direction  (HCYL); 


Vertical  Elliptical  Cylinder 
Oriented  Perpendicular  To 
Flow  With  L>  W  Or  W>L 
(VCYL); 


Block  (BLOCK): 


Streamline  Half  Bodies 
(FAIRING); 


Launch  Shoe  (SHOE); 

1)  Base  Member 

2)  Vortical  Member 

3)  Top  Member 


V 


TOP  VIEW 


3-D  VIEW 


CO 

1 

_ ^ 

i 

Launch  Lug  (LUG): 

1)  Base  Member 

2)  &  3)  Vertical  Members  ^ 

4)  Top  Member 

L-Length 


W-Wldth  H-Helght 


Note;  Length,  width,  height,  and  offset  must  be  Input  lor  each  member  of  launch  lug  and 
launch  shoe  types 

Note;  Offset  is  the  perpendicular  distance  from  the  missile  mold  line  to  the  bottom  of  the 
protuberance  or  protuberance  member 

Figure  8  Protuberance  Shapes  Available 
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CASBID  PROTUBBRAMCB  EXAMPLE  CASE 
DIM  IN 
NO  LAT 

$RBFQ  XCG=39.0,$ 

$FLTCON  NMACH=3 . ,  MACH=0 . 4, 0 . 8, 2 . 0, 

REN=3 . B06, 3 .B06, 3 .B06, ALT=0 . 0, 

NALPHAsS  .  ,  ALPHA;=-8  .,-4.,0.,4.,8.,$ 
$AXZBOD  TNOSB=OGIVB, LNOSX=12 . 0,DNOSE=12 .0, 
LCXMTtlsS4 . 0,DCEM7R^12 . 0, 

TAFT«CONE ,  IiAl‘T*12 . 0 ,  DAFT=6 . 0 ,  DEXIT=5 . 0 ,  $ 
$PROTUB  NPROT»4., 

PTYPE^rAIRING,  VCYX, SHOE, BLOCK, 

XPROT“14 . , 22 . ,  39 . , 56 . , 

NL0C»2 . , 4 . , 2 . , 1 . , 

LPROT=5 . ,  1 . ,  10 . ,  10 . ,  10 . ,  0 . 5, 

WPROT=2  ,,1,,4.,0.25,1.,1., 

HPROT=2 .,0.5,0.1,0.75,0.25,0.25, 
OPROT=0.,0.,0., 0.1,  0.85,0. ,$ 

$FINSET1  SSPANsO.O, 9.0,CHORD=14.0, 8.0, 

XLES64 . 0 , SNEBPsQ . 0, 8TA=1 . 0 , NPANELs4 . , 
PHir=45  . ,  13.5 . ,  225 . ,  315 . ,  $ 

PRINT  GEOM  BODY 
PRINT  AERO  BODY 
SAVE 

NEXT  CASE 


NOTE;  Length,  Width,  and  Haight  it  Input  tor  aach  mambar  of  tha  launch  ahoa 

Figure  9  Protuberance  F.xample  Input  File 
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NAMELIST  FINSETn 

NOMINAL  INPUTS 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

DERATION 

UNITS 

DEFAULT 

SECTYP 

TYPE  OF  SECTION  TO  BE  DEFINED;  HEX. 

- 

HEX 

sspanO 

NACA.  ARC.  OR  USER 

10 

SEMI-SPAN 

L 

- 

CHORD 

10 

PANEL  CHORD  LENGTH  AT  EACH  SSPAN 

L 

- 

XLE 

10 

DISTANCE  FROM  NOSE  TIP  TO  CHORD 

L 

0.0 

SWEEP© 

LEADING  EDGE  AT  EACH  SSPAN 

10 

SWEEPBACK  ANGLE  AT  EACH  SSPAN 

DEG 

0.0 

STA 

10 

CHORD  STATION  USED  IN  MEASURING 

- 

1.0 

LER© 

SWEEP  AT  EACH  SSPAN  (O  O^^LEADING 
EDGE.  I  .OrTRAILING  EDGE) 

10 

PANEL  LEADING  EDGE  RADIUS  AT  EACH 

L 

• 

SPAN  STATION 

NPANEL 

• 

NUMBER  OF  PANELS  IN  SET  (1  -8) 

• 

4 

PHIF 

8 

ROLL  ANGLE  OF  EACH  FIN  MEASURED 

DEG 

@ 

CLOCKWISE  FROM  TOP  VERTICAL 
CENTER  LOOKING  FORWARD 

QAM 

3 

DIHEDRAL  OF  EACH  FIN  (POSITIVE  WHEN 

DEG 

0.0 

PHIF  IS  INCREASED.  SEE  FIG.  12) 

SKEW 

- 

ANGLE  BETWEEN  THE  Y  AXIS  AND  THE 

DEG 

0.0 

FIN  HINGE  LINE  (POSH  IVE  SWEPT  BACK) 

(C  IF  SSPAN(1)-0.0.  INPUTS  ARE  RELATIVE  TO  ROOT  CHORD  NOT  BODY  CENTERLINE 
®  IF  USING  SWEEP.  SPECIFY  ONLY  XLE{1);  IF  USING  XLE  DO  NOT  SPECIFY  SWEEP 
(3)  NOT  REQUIRED  FOR  NACA  AIRFOILS.  REQUIRED  FOR  USER  AIRFOILS 
®  IF  PHIF  NOT  INPUT  THE  NUMBER  OF  PANELS  ARE  EVENLY  SPACED  ABOUT  THE  BODY, 


Figure  10a  Fin  Geometry  Inputs  -  Nominal 


NAMELIST  FINSETn 

OPTIONAL  INPUTS 


ARRAY 

DIMENSION 

DERNITION 

UNITS 

DEFAULT 

ZUPPER 

10 

MAXIMUM  THICKNESS  TO  CHORD  RATIO 

OF  UPPER  SURFACE 

- 

0.025 

ZLOWER 

10 

MAXIMUM  THICKNESS  TO  CHORD  RATIO 

OF  LOWER  SURFACE 

ZUPPER 

LMAXU 

10 

FRACTION  OF  CHORD  FROM  SECTION 
LEADING  EDGE  TO  MAXIhtUM 

THICKNESS  OF  UPPER  SURFACE 

0.5 

LMAXL 

10 

FRACTION  OF  CHORD  FROM  SECTION 
LEADING  EDGE  TO  MAXIMUM 

THICKNESS  OF  LOWER  SURFACE 

LMAXU 

LFLATU 

10 

FRACTION  OF  CHORD  OF  CONSTANT 
THICKNESS  SECTION  ON  UPPER 
SURFACE 

0.0 

LFLATL 

10 

1 

FRACTION  OF  CHORD  OF  CONSTANT 
THICKNESS  SECTION  ON  LOWER 
SURFACE 

LFLATU 

NOTE;  THESE  PARAMETERS  MUST  BE  INPUT  AT  EACH  SPAN  STATION 


Figure  10b  Fin  Geometry  Inputs  -  Optional 
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NAMELIST  FINSETn 

INPUTS  FOR  "USER"  SECTIONS 


VARIABLE 

ARRAY 

DERNITION 

NAME 

DIMENSION 

UNITS 

DEFAULT 

XCORD 

50 

CHORD  STATION,  FRACTION  OF  CHORD 

MEAN  0 

FROM  LEADING  EDGE 

50 

DISTANCE  BETWEEN  THE  MEAN  LINE  AND 

- 

. 

CHORD  LINE  AT  EACH  XCORD. 

FRACTION  OF  CHORD 

THICK  0 

50 

THICKNESS  TO  CHORD  RATIO  AT  EACH 

XCORD 

YUPPER  0 

50 

UPPER  SURFACE  COORDINATES. 

. 

- 

FRACTION  OF  CHORD.  AT  EACH  XCORD 

YLOWER  (T 

)  50 

LOWER  SURFACE  COORDINATES. 

- 

1 

FRACTION  OF  CHORD.  AT  EACH  XCORD 

NOTE:  ALL  VARIABLES  ARE  EXPRESSED  AS  FRACTIONS  OF  CHORD 
LEADING  EDGE  RADIUS  (VARIABLE  LER)  MUST  BE  DEFINED 

®  EITHER  MEAN  AND  THICK  OR  YUPPER  AND  YLOWER  ARE  REQUIRED 


+Y 


Figure  10c  Fin  Geometry  Inputs  -  User  Airfoils 
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MISSILE  BODY  CENTERUNE 


0 


NOTE;  XLE  IS  MEASURED  FROM  NOSE  TIP 

IF  SSPAN<1)  IS  INPUT  AS  ZERO,  SSPAN  INPUTS 
ARE  RELATIVE  TO  BODY  SURFACE  MOLD  LINE 

Figure  11  Selecting  Panel  Break  Points 
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t‘ 

TOP  VERTICAL 
CENTER 


PANEL  2 


VIEW  LOOKINQ  FORWARD 
FROM  BEHIND  MISSILE 


Figure  12  Fin  Numbering  and  Orientation 
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VIEW  LOOKING  FORWARD 
FROM  BEHmO  MISSILE 


PHI  IS  THE  BODY  ROLL  ANGLE 
PHIF  IS  THE  HN  PANEL  ROLL  ANGLE 


Figure  13  Roll  Attitude  vs  Fin  Orientation 


NAMELIST  DEFLCT 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

OEFINI-nON 

UNITS 

DEFAULT 

DELTA1 0 

8 

DEFLECTION  ANGLES  FOR  EACH  PANEL  IN 
FIN  SET  1  (SUBSCRIPT  IS  FIN  NUMBER) 

DEO 

0. 

DELTA20 

8 

DEFLECTION  ANGLES  FOR  EACH  PANEL  IN 
FIN  SET  2  (SUBSCRIPT  IS  FIN  NUMBER) 

DEG 

0. 

DELTA30 

8 

DEFLECTION  ANGLES  FOR  EACH  PANEL  IN 
FIN  SET  3  (SUBSCRIPT  IS  FIN  NUMBER) 

DEG 

0. 

8 

DEFLECTION  ANGLES  FOR  EACH  PANEL  IN 
FIN  SET  4  (SUBSCRIPT  IS  FIN  NUMBER) 

DEG 

0. 

XHINGE 

4 

DISTANCE  FROM  ROOT  CHORD  LEADING 
EDGE  TO  PANEL  HINGE  LINE  FOR  EACH 
FIN  SET 

L 

CR/2  0 

SKEW 

4 

SWEEPBACK  OF  HINGE  LINE  FOR  EACH 

FIN  SET 

DEG 

0. 

<D  PANEL  NUMBERING  IS  SHOWN  IN  FIGURE  12 
(S  DEFAULT  IS  AT  ONE-HALF  THE  EXPOSED  ROOT  CHORD 

NOTE;  A  POSITIVE  DEFLECTION  ANGLE  PRODUCES  A  NEGATIVE  BODY  AXIS 
ROLLING  MOMENT  AT  ZERO  ANGLE  OF  ATTACK 


Figure  14  Panel  Deflection  Inputs 
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NAMELIST  TRIM 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

DERATION 

UNITS 

DEFAULT 

SET 

FIN  SET  TO  BE  USED  FOR  TRIM 

1. 

PANL1 ' 

.TRUE.  IF  PANEL  TO  BE  USED 

.FALSE. 

PANL2 

.TRUE.  IF  PANEL  TO  BE  USED 

•FALSE. 

PANL3 

•0 

.TRUE.  IF  PANEL  TO  BE  USED 

.FALSE. 

PANL4 

.TRUE.  IF  PANEL  TO  BE  USED 

.FALSE. 

PANL5 

•TRUE.  IF  PANEL  TO  BE  USED 

FALSE. 

PANL6 

.TRUE.  IF  PANEL  TO  BE  USED 

.FALSE. 

PANL7 

.TRUE.  IF  PANEL  TO  BE  USED 

.FALSE, 

PANL8, 

.TRUE.  IF  PANEL  TO  BE  USED 

.FALSE. 

DELMIN 

MINIMUM  NEGATIVE  DEFLECTION 

DEG 

-25. 

DELMAXrM 

MAXIMUM  POSITIVE  DEFLECTION 

DEG 

+20. 

ASYM 

4 

.TRUE.  IF  PANEL  WITH  SUBSCRIPT  IS  TO 

BE  DEFLECTED  OPPOSITE  TO  NORMAL 
SIGN  CONVENTION  (ASYMMETRIC 
DEFLECTIONS) 

FALSE. 

(P  DEFAULTS  APPLY  ONLY  IF  ALL  PANL#  DATA  ARE  NOT  INPUT  OR  .FALSE. 
®  BOTH  DELMIN  AND  DELMAX  MUST  BE  SPECIFIED 


Figure  15  Trim  Inputs 
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NAMELIST  INLET 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

DERNITION 

UNITS 

DEFAULT 

ININ 

- 

NUMBER  OF  INLETS  (MAXIMUM«20) 

- 

- 

INTYPE 

- 

TYPE  OF  INLET:  2DTOP,  2DSIDE,  OR  AXI  0 

- 

- 

XINLT 

■ 

LONGITUDINAL  DISTANCE  FROM  NOSE  TIP 
TO  INLET  LEADING  EDGE 

L 

• 

XDIV 

LONGITUDINAL  DISTANCE  FROM  INLET 
LEADING  EDGE  TO  DIVERTER 

LEADING  EDGE 

L 

HDIV 

- 

HEIGHT  OF  DIVERTER  LEADING  EDGE 

L 

- 

LDIV 

- 

LENGTH  OF  DIVERTER 

L 

- 

|PHI® 

20 

INLET  ROLL  ORIENTATIONS 

DEG 

- 

X  ' 

5 

INLET  LONGITUDINAL  POSITIONS  RELATIVE 
TO  INLET  LEADING  EDGE 

L 

• 

H® 

5 

INLET  HEIGHTS  AT  THE  LONGITUDINAL 
POSITIONS 

L 

• 

W©J 

5 

INLET  WIDTHS  AT  THE  LONGITUDINAL 
POSITIONS 

L 

* 

■cover 

• 

IF  COVER-.TRUE.  INLETS  ARE  COVERED 

• 

FALSE. 

RAMP 

• 

EXTERNAL  COMPRESSXDN  RAMP  ANGLE 

DEO 

- 

ADD 

• 

IF  ADD-.TRUE.  INLET  ADDITIVE  DRAG  IS 
CALCULATED 

* 

FALSE. 

MFR0 

20 

MASS  FLOW  RATIO  FOR  EACH  MACH 
NUMBER  (OsMFRi  1.0) 

■ 

(D  2DTOP:  TWO  DIMENSIONAL  TOP  MOUNTED,  2DSIDE:  TWO  DIMENSIONAL  SIDE  MOUNTED, 
AXI:  AXISYMMETRIC 

<J)  ROLL  POSITIONS  FROM  TOP  VERTICAL  CENTER.  SAME  CONVENTION  AS  FIN  ROLL 
POSITIONS. 

<3)  SPECIFY  X,  H,  AND  W  AT  FIVE  INLET  LOCATIONS;  1)  LEADING  EDGE,  2)  COWL  LIP  LEADING 
EDGE,  3)  MIDBODY  START,  4)  START  OF  BOATTAIL,  5)  END  OF  BOATTAIL 
®  NOT  REOUIRED  IF  INTYPE-AXI 
®  IF  INTYPE.AXI,  W-DIAMETER 

(B>  SPECIFY  MASS  FLOW  RATIO,  MFR,  FOR  EACH  FREE  STREAM  MACH  NUMBER  OR  VELOCITY 
GIVEN  IN  SFLTCON  (ONLY  REQUIRED  IF  ADD  =  TRUE.) 


Figure  16  Inlet  Geometry  Inputs 
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SIDE  VIEW 


(1) 

X(1)-0 

H(1)-0 


BOTTOM  VIEW 


TOP  VERTICAL 
CEMTER 

♦ 

O- 


INPUT  INLET  LONOrrUOINAL  COORDINATE  (X),  HEIGHT  (H), 
AND  WIDIH  (W)  AT  EACH  STATION 

STATION  INLET  POSITION 

(t)  LEADING  EDGE  (X(1)&H(1)..0) 

(2)  COWL  LIP  LEADING  EDGE 

(3)  MIDBODY  START 

(4)  BOATTAIL  START 

(5)  BOATTAIL  END 


INLET  ALONE 


FRONT  VIEW 


-P 


1“^  I  I 

(1)  (2)  (3) 


NOTES; 

•  INLET  ROLL  ORIENTATION  IS  SAME  CONVENTION  AS  FIN  ROLL  ORIENTATION. 

•  RAMP  IS  THE  EXTERNAL  COMPRESSION  RAMP  ANGLE  (SHOWN  SHADED  IN  THE  SIDE  VIEW) 

•  HEIGHT  OF  THE  DIVERTER  IS  SPECIFIED  AT  THE  DIVERTER  LEADING  EDGE 

•  THE  DIVERTER  WIDTH  IS  EQUAL  TO  THE  INLET  WIDTH  AT  LDIV 

■  IF  INLET  IS  COVERED  (COVER-  TRUE.)  A  PLUG  IS  PLACED  BETWEEN  STATIONS  1  AND  2 
FLUSH  WITH  THE  INLET  FACE 


Figiire  17a  Top-Mounted  2-D  Inlet /Diverter  Geometry 
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TOP  VIEW 


H(1)-0 


SO" - 


TOP  VERTICAL 
CEHTER 


INPUT  INLET  LONOITUDINAL  COORDINATE  (X),  HEIGHT  (H), 
AND  WIDTH  (W)  AT  EACH  STATION 

STATION  INLET  POSITION 

(1)  LEADING  EDGE  (X(1)&H(1)-0) 

(2)  COWL  LIP  LEADING  EDGE 

(3)  MIDBODY  START 

(4)  BOATTAIL  START 

(5)  BOATTAIL  END 


180" 


INLET  ALONE 


FRONT  VIEW 


NOTES; 

.  INLET  ROLL  ORIENTATION  IS  SAME  CONVENTION  AS  FIN  ROU  ORIENTATION 

•  RAMP  IS  THE  EXTERNAL  COMPRESSION  RAMP  ANGLE  (SHOWN  SHADED  IN  THE  SIDE  VIEW) 
.  HEIGHT  OF  IHfc  DIVERTER  IS  SPECIFIED  AT  THE  DIVERTER  LEADING  EDGE 

.  THE  DIVERTER  WIDTH  IS  EQUAL  TO  THE  INLET  WIDTH  AT  LDIV 

•  IF  INLET  IS  COVERED  (COVER-TRUE  )  A  PLUG  IS  PLACED  BETWEEN  STATIONS  1  AND  2 
FLUSH  WITH  THE  INLET  FACE 


Figure  17b  Side-Mounted  Z-D  Inlet/ Diverter  Geometry 


56 


«UnF  UIFUU 


(1) 

X(1).0 

W(1)-0 


BOTTOM  VIEW 


I  I  I 

(1)  (2)  (3) 


TOP  VERTICAL 
CEKTER 

♦ 

^0* 


INPUT  INLET  LONGITUDINAL  COORDINATE  (X),  AND 
DIAMETER  (W)  AT  EACH  STATION 


/  ✓ 

PKI(l>.t20i 


STATION 

(1) 

(2) 

(3) 

(4) 

(5) 


INLET  POSITION 

TIP  (X(1)&  W(1)-0) 
COWL  LIP  LEADING  EDGE 
MIDBODY  START 
BOATTAIL  START 
BOATTAIL  END 


PHI<2).2S0.0' 


INLET  ALONE 


|Z  inlets  rotated  at  120*  AND  iftn 


FRONT  VIEW 


inlet  FACE 


"Kl  I 

(1)  (2)  (3) 


NOTES; 

.  INLET  ROLL  ORIENTATION  IS  SAME  CONVENTION  AS  FIN  ROLL  ORIENTATION 

•  RAMP  IS  THE  EXTERNAL  COMPRESSION  CONE  HALF-ANGLE  (SHOWN  SHADED 
IN  THE  SIDE  VIEW) 

•  HEIGHT  OF  THE  DIVERTER  IS  SPECIFIED  AT  THE  DIVERTER  LEADING  EDGE 

•  THE  DIVERTER  WIDTH  IS  EQUAL  TO  THE  INLET  DIAMETER  AT  LDIV 

.  IF  INLET  IS  COVERED  (COVER-.TRUE  )  A  PLUG  IS  PLACED  BETWEEN  STATIONS  1  AND  2 
FLUSH  WITH  THE  INLET  FACE 


Figure  17c  Axisymmetric  Inlet /Diverter  Geonretry 


J 
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NAMELIST  EXPR 


VARIABLE 

NAME 

ARRAY 

DIMENSION 

OERNITION 

UNITS 

DEFAULT 

MACH 

- 

MACH  NUMBER 

- 

- 

NALPHA 

- 

NUMBER  OF  ANGLES  OF  ATTACK  (2-20) 

- 

- 

ALPHA 

20 

ANGLES  OF  ATTACK  FOR  DATA 

DEG 

- 

SREF 

- 

REFERENCE  AREA  FOR  DATA 

L'L 

LREF 

• 

LONGITUDINAL  REFERENCE  LENGTH 

FOR  DATA 

L 

LATREF 

LATERAL  REFERENCE  LENGTH  FOR  DATA 

L 

LREF 

XCQ 

- 

LONGITUDINAL  C.Q.  FOR  DATA 

L 

0. 

ZCG 

VERTICAL  C.G.  FOR  DATA 

L 

0. 

CONF 

CONFIGURATION  FOR  DATA 

SELECT  ONE  OF  THE  FOLLOWING 

BODY  -  BODY 

F1  -  WING 

F2  -  TAIL 

F3  -  THIRD  FIN  SET 

F4  -  FOURTH  FIN  SET 

8F1  -  BODY-WING 

BFI2- BODY-2  FIN  SETS 

BF123- BODY-3  FIN  SETS 

BF1 234  -  BOOY-4  FIN  SETS 

CN 

20 

Cn  data  VS  ALPHA 

• 

CM 

20 

Cm  data  VS  ALPHA 

- 

CA 

20 

C*  DATA  VS  ALPHA 

- 

CY 

20 

Cr  DATA  VS  ALPHA 

« 

CSN 

20 

Cn  DATA  VS  ALPHA 

- 

CSL 

-  20  _. 

Cl  DATA  VS  ALPHA 

- 

(D  DEFAULT  IS  BODY  MAXIMUM  CROSS-SECTIONAL  AREA.  IF  NO  BODY  IS  INPUT. 
MAXIMUM  FIN  PANEL  AREA  IS  USED. 

®  DEFAULT  IS  BODY  MAXIMUM  DIAMETER.  IF  NO  BODY  IS  INPUT.  MAXIMUM  FIN  PANEL 
MEAN  GEOMETRIC  CHORD  IS  USED. 


Figure  18  Experimental  Data  Inputs 
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Table  5  Common  Block  DUMP  and  WRITE  Names 


COMMON  BLOCK 


DUMP  NAME 


WRITE  NAME 


ABOOIN 

BDWORK 

CASEIO 

CON9V 

DBODY 

DB1 

0B12 

0B123 

DB1234 


OFLAGS 

DUMPF 


FSET1 
FSET2 
FSET3 
FSET4 
FI  WORK 
F2WORK 
P3WORK 
F4WORK 
GEOBOD 
QEOFS1 
GEOFS2 
QEOFS3 
GEOFS4 
INCID 
INLETN 
INLTO 
INPCON 
LOGIC 
PAERO 
REFQN 
SBODY 
SB1 
SB12 
SB123 
SB1234 
SFIN1 
SFIN2 
SFIN3 
SFIN4 
THERY 
TOTALC 
TRACE 
TRIMD 
TRIMIN 
UTRIMD 


BDIN 

BDWK 


DBOD 

DB1 

DB12 

0B13 

DB14 


FLT 

F1IN 

F2IN 

F3IN 

F4IN 

F1WK 

K2WK 

F3WK 

F4WK 

GEOB 

F1GM 

F2GM 

F3GM 

F4GM 


REFQ 

SBOD 

SB1 

SB12 

SB13 

SB14 

SF1 

SF2 

SF3 

SF4 


ABODIN  or  EBODIN 
BDWORK 
CASEID 
CONST 
OBODY 
DB1 
DB12 
DB123 
DB1234 
DESIG 
DFIN1 
DFIN2 
0FIN3 
DFIN4 
DFLAGS 
DUMPF 
FLC 
FSET1 
FSET2 
FSET3 
FSET4 
FI  WORK 
F2WORK 
F3WORK 
F4WORK 
GEOBOD 
GEOFS1 
GEOFS2 
GEOFS3 
GEOFS4 
INCID 
INLETN 
INLTD 
INPCON 
LOGIC 
PAERO 
REFQN 
8BODY 
SB1 
SB12 
SD123 
SB1234 
SFIN1 
SFIN2 
SFIN3 
SFIN4 
THERY 
TOTALC 
TRACE 
TRIMD 
TRIMIN 
UTRIMD 
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Table  6  Airfoil  Designation  Using  the  NACA  Control  Card 


INPUT  NACA 
DESIGNATION 

NACA  SERIES 
AIRFOIL 

RESTRICTIONS 

0012.25 

4-Olgit 

None.  Fractional  thickrtess  may  be 
specified. 

23118.50 

5-Digit 

None.  Fractional  thickness  may  be 
specified. 

2406-32 

4-DigH 

modified 

Sixth  digit  specifies  position  of 
maximum  thickness,  (%chord/10),  and 
must  be  a  2, 3, 4, 5,  or  6. 

43006-65 

5-Oigit 

modified 

Seventh  digit  specifies  posHlon  of 
maximum  thickness,  (%chord/10),  and 
must  be  a  2, 3, 4, 5,  or  6. 

16-212.25 

1 -Series 

Second  digit  specifies  location  of 
minimum  pressure,  (%chord/10),  and 
must  be  a  6, 8,  or  9.  Fractional 
thickness  may  be  specified. 

64-005 

64-205  AsO.6 
63A005 

652A215  AsO.8 
65,2A215  A=0.8 

6-Serle8 

Second  digit  specifies  location  of 
minimum  pressure,  (%chord/10),  and 
must  be  a  3, 4, 5,  or  6.  The  mean  line 
parameter  (Asxx)  must  be  a  decimal 
between  0.1  and  1.0  (Default  isl.O). 

See  Note  1. 

3-30.0-2.5-40.1 

A  B  C  0 

Supersonic 

See  Note  2. 

A  -  Section  type;  1=Double  Wedge 
2sCircular  Arc 
3sHexagonal 

B  -  Distance  from  leading  edge  to 
position  of  maximum  thickness,  %  of 
chord. 

C  -  Maximum  thickness,  %  of  chord, 

D  -  For  hexagonal  sections,  length  of 
surface  of  constant  thickness,  %  of 
chord. 

Note  1 .  The  program  does  not  distinguish  between  a  64,2-220  and  a  64-220 
specification.  The  difference  In  coordinates  between  the  two  is  negligible. 


Note  2.  All  parameters  can  be  expressed  to  0.1%.  The  delimiter must  be  used. 
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Figure  20  "Configuration  Incrementing”  Case  Set-up 
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4.0  OUTPUT  DESCRIPTION 


This  section  describes  the  types  of  output  available  from  the  code.  In 
most  cases  the  available  output  is  user  selectable,  that  is,  it  is  not  normally 
provided  and  must  be  specifically  requested  using  a  specialized  control  card. 
This  feature  permits  the  user  to  tailor  the  code  output  to  fit  his  particular 
application  without  extensive  reprogramming.  This  allows  him  to  find  the 
output  that  he  is  interested  in  without  having  to  wade  through  output  that 
does  not  interest  him. 

The  following  four  types  of  output  are  available  from  the  code: 

•  Nominal  output  -  This  output  is  always  provided  by  the  code 
and  consists  of  output  from  the  input  error  checking  module 
(CONERR),  a  listing  of  the  inputs  for  each  case,  and  the  final 
aerodynamic  results  for  the  configuration. 

•  Partial  output  -  This  output  details  the  configuration 
geometry  and  the  intermediate  aerodynamic  calculations. 

Special  control  cards  are  available  so  that  the  user  can  select 
the  quantity  and  types  of  output  desired. 

•  External  data  files  -  This  output  permits  the  user  to  create 
external  data  files  which  can  be  used  in  post-processing 
programs,  such  as  plotting  or  trajectory  programs.  Both  fixed 
and  user  defined  format  data  files  can  be  created  with  the 
addition  of  simple  control  cards. 

•  Array  dumps  and  extrapolation  messages  -  This  output 
permits  the  user  to  print  internal  data  arrays  (DUMP)  or  to 
investigate  design  chart  extrapolations  during  the  course  of 
the  case  execution  (PRINT  EXTRAP). 

The  remainder  of  the  section  describes  each  of  these  output  data. 
Examples  of  each  output  page  are  also  included  and  were  created  from  the 
example  problems,  described  in  the  Appendix,  which  can  be  used  as  a  model 
for  setting  up  another,  similar  configuration  or  be  used  as  a  means  to  check 
the  proper  operation  of  the  code. 

4.1  Nominal  Output 

Without  the  use  of  any  program  options  the  code  will  provide  three 
types  of  output.  First,  an  analysis  by  the  input  error  checking  routine  is 
provided.  It  lists  all  input  cards  provided  by  the  user  and  identifies  any  input 
errors  detected.  Second,  a  listing  of  all  input  cards,  grouped  by  case,  are 
provided;  included  in  this  output  is  an  error  analysis  from  the  major  input 
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error  routine  MAJERR.  Finally,  the  total  configuration  aerodynamics  are 
provided  in  summary  form;  one  page  of  aerodynamic  output  is  supplied  for 
each  Mach  number  specified.  The  MAJERR  results  and  the  total  configuration 
aerodynamics  results  are  listed  in  succession  for  each  case. 

4.1.1  Input  Error  Checking 

The  purpose  of  the  input  error  checking  module  is  to  provide  single 
pass  error  dtecking  of  all  inputs.  If  an  error  is  detected,  it  is  identified  and  an 
appropriate  error  message  provided.  The  error  messages  are  designed  to  be 
self-explanatory.  In  some  cases,  errors  are  automatically  corrected  by  the 
routine,  although  the  routine  was  not  designed  to  be  a  comprehensive  error 
correction  utility. 

The  following  errors  are  automatically  corrected  by  the  code: 

•  No  terminating  comma  on  a  namelist  input  card 

•  No  terminating  "$”  or  "$END"  on  a  namelist  input  (”&"  on 
IBM  systems) 

•  No  terminating  NEXT  CASE  for  the  case  inputs  for  single 
case  or  last  case  inputs. 

Errors  detected  by  the  error  checking  routine  are  considered  either 
"FATAL"  or  'NON-FATAL",  a  "FATAL"  error  is  one  which  will  cause  the 
code  to  terminate  execution  abnormally;  examples  of  "FATAL"  errors  include 
incorrect  spelling  of  any  namelist  name,  incorrect  spelling  of  any  variable 
name,  and  any  drastic  input  error  in  a  namelist  input,  such  as  leaving  out  an 
equals  sign  in  a  constant  definition.  All  "FATAL"  errors  are  clearly  identified 
on  the  output.  A  "NON-FATAL"  error  is  one  which  will  not  cause  the 
program  to  terminate  execution;  an  example  of  a  "NON-FATAL"  error  is 
leaving  off  the  decimal  point  on  numeric  constants  all  Missile  Datcom  inputs 
are  either  REAL  or  LOGICAL  regardless  of  the  variable  name  assigned. 
"NON-FATAL"  errors  will  not  cause  the  code  to  stop  execution,  whereas, 
"FATAL"  errors  will  cause  the  code  to  stop  execution  after  input  error 
checking  has  been  completed. 

An  example  output  from  CONERR  is  shown  in  Figure  21.  This  figure 
illustrates  the  array  of  input  errors  checked  by  CONERR.  Several  additional 
features  of  the  output  are  as  follows: 

•  All  user  defined  input  cards  are  assigned  a  sequential  "line 
number".  This  serves  to  identify  user  inputs  from  the  code 
generated  inputs  (all  code-created  input  cards  are  not 
identified  with  a  "line  number").  This  scheme  also  permits 
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the  user  to  quickly  identify  those  input  cards  in  error  so  that 
efficient  correction  of  input  errors  can  be  performed. 

•  All  input  cards  are  listed  as  input  by  the  user.  To  the  right  of 
each  input  card  is  a  listing  of  any  errors  encountered  in 
processing  that  card.  If  no  such  error  message  appears  then 
the  input  was  interpreted  as  being  correct. 

•  In  many  cases  alphanumeric  constants  are  available  (see 
Table  3).  Hence  the  user  does  not  need  to  memorize  a 
numeric  scheme  of  "flags".  Since  some  computers  do  not 
recognize  alphanumeric  constants  as  namelist  constants,  they 
are  automatically  converted  by  the  code  to  their  numeric 
equivalent.  A  message  is  printed  to  identify  the  substitutions 
performed.  The  example  input  in  Figure  21  shows 
replacements  for  CONE  and  OGIVE. 

In  order  to  permit  colurrm  independent  inputs  the  code  will 
automatically  adjust  some  of  the  input  cards  to  begin  in  columns  1  or  2.  All 
control  cards  will  be  automatically  shifted  to  start  in  column  1;  all  namelists 
which  begin  in  column  1  will  be  shifted  to  column  2.  If  any  input  card  cannot 
be  shifted  to  conform  to  this  scheme,  an  error  message  will  be  produced.  As  a 
general  rule,  column  80  of  namelist  inputs  should  be  left  blank  so  that  the 
code  can  shift  the  card  image,  if  necessary. 

4.1.2  Listing  of  Case  Input  Data 

Figure  22  shows  the  first  page  of  outputs  for  a  case  without  CONERR 
detected  errors.  Then  Figure  23  shows  the  next  page  of  output  which  lists  all 
input  cards  for  the  case  (down  to  the  NEXT  CASE  control  card).  If  the  input 
for  a  case  is  from  a  previous  case  (through  use  of  the  SAVE  control  card)  only 
the  new  case  inputs  are  listed.  All  saved  inputs  are  not  repeated  in  subsequent 
case  input  summaries. 

After  the  case  data  have  been  read,  the  data  setup  for  the  case  is 
analyzed  by  the  case  major  error  checking  module  (MAJERR).  The  purpose  of 
this  second  error  checking  is  to  insure  that  the  data  input,  although  syntax 
error  free,  properly  defines  a  case  to  be  run.  Examples  of  errors  detected  in 
MAJERR  include  valid  flight  condition  inputs,  valid  reference  condition 
inputs,  and  that  geometry  has  been  defined.  In  most  cases  errors  detected  by 
MAJERR  are  corrected  with  assumed  defaults.  If  any  MAJERR  error  message 
is  produced,  the  user  should  verify  the  "fix-up"  taken  by  the  code.  In  some 
cases  a  "nx-uf>  is  not  possiOle;  an  appropriate  error  message  and  a  suggestion 
for  correcting  the  error  is  provided.  If  a  "fix-up"  is  not  possible  the  case  will 
not  run. 
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4.1.3  Case  Total  Configuration  Aerodynamic  Output  Summary 


As  shown  in  Figure  24,  the  total  configuration  aerodynamics  are 
provided  in  compact  form  for  easy  review.  The  aerodynamics  are 
summarized  as  a  function  of  angle  of  attack  (ALPHA)  in  the  user  specified 
system  of  units,  the  nomenclature  is  as  follows: 

CN  -  Normal  force  coefficient 

CM  -  Pitching  moment  coefficient 

CA  -  Axial  force  coefficient 

CY  -  Side  force  coefficient 

CLN  -  Yawing  moment  coefficient 

CLL  -  Rolling  moment  coefficient 

CNA  -  Normal  force  coefficient  derivative  with  ALPHA 

CMA  -  Pitching  moment  coefficient  derivative  with 

ALPHA 

CYB  -  Side  force  coefficient  derivative  with  BETA 

CLNB  -  Yawing  moment  coefficient  derivative  with  BETA 

CLLB  -  Rolling  moment  coefficient  derivative  with  BETA 

CL  -  Lift  coefficient 

CD  -  Drag  coefficient 

CL/CD  -  Lift  to  drag  ratio 

XCP  -  Center  of  pressure  from  the  moment  reference 

center  divided  by  reference  length 

All  coefficients  are  based  upon  the  reference  areas  and  lengths  specified 
at  the  top  of  the  output  page.  The  derivatives  CNA  and  CMA  are  computed 
by  numeric  differentiation  of  the  CN  and  CM  curves,  respectively;  precise 
derivatives  are  only  obtained  when  the  angle  of  attack  range  specified  is 
narrow.  The  derivatives  CYB,  CLNB  and  CLLB  are  determined  by  perturbing 
the  sideslip  angle  by  one  degree,  recalculating  the  configuration  forces  and 
moments,  and  then  differencing  with  the  user  specified  orientation.  Hence, 
the  longitudinal  and  lateral  derivatives  will  probably  not  be  numerically 
identical  for  thosr  conditions  which  should  produce  identical  results  if  they 
were  both  calculated  by  the  same  method. 

A  significant  decrease  in  computational  time  is  realized  when  the 
calculation  of  lateral-directional  derivatives  are  suppressed  using  the  control 
card  NO  LAT.  For  these  cases,  the  CYB,  CLNB,  and  CLLB  data  fields  are  filled 
with  blanks. 

When  selecting  TRIM,  the  output  is  provided  in  a  form  similar  to 
Figure  25.  When  running  a  trii  t  case  the  derivatives  due  to  ALPHA  and 
BETA  ;  .e  not  available.  The  panels  which  were  deflected  to  trim  the 
configuration  are  indicated  by  the  "VARIED"  citation  next  to  them. 
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The  format  for  the  values  of  the  numbers  in  the  printed  output  has 
been  assumed  based  on  typical  magnitudes  for  missile  aerodynamic 
coefficients.  In  some  cases,  a  user  specified  reference  area  and/or  length  will 
cause  the  results  to  underflow  or  overflow  the  format  selected.  For  these  cases 
the  user  should  adjust  his  reference  quantities  by  powers  of  ten  to  get  the  data 
to  fit  the  format  specified. 

4.2  Partial  Output 

Partial  output  consists  of  geometry  calculation  details,  intermediate 
aerodynamic  results,  or  auxiliary  data,  such  as  pressure  distributions.  Each  of 
these  output  types  are  printed  through  the  addition  of  control  cards  input  for 
each  case.  In  all  cases,  partial  output  requested  for  one  case  is  not 
automatically  selected  for  subsequent  cases,  and  the  control  cards  must  be  re¬ 
input.  This  permits  the  user  to  be  selective  on  the  amount  and  types  of 
output  desired. 

A  special  control  card  PART  permits  the  user  to  request  all  geometric 
and  aerodynamic  partial  output.  Due  to  the  amount  of  output  produced,  this 
option  should  be  used  sparingly  or  when  details  of  the  calculations  are 
desired. 

There  is  one  geometry  for  which  no  partial  output  is  obtained.  This  is 
when  bodies  with  arbitrary  cross  sections  are  input. 

The  following  paragraphs  describe  the  output  received  when  partial 
output  is  requested. 

4.2.1  Geometric  Partial  Output 

Details  of  the  geometry  are  provided  when  the  PART  or  PRINT  GEOM 
control  cards  are  included  in  the  case  inputs.  Figure  26  shows  the  output 
created  when  the  PRINT  GEOM  BODY  control  card  is  used.  Detailed  are  the 
results  of  the  geometric  calculations  for  the  body.  Included  are  such  items  as 
planform  area,  surface  (wetted)  area,  and  the  mold  line  contour. 

If  fins  are  present  on  the  configuration,  tv^o  types  of  fin  geometry  data 
are  produced  when  PRINT  GEOM  FINl  or  PART  is  requested.  As  shown  in 
Figure  27,  the  description  of  the  panel  airfoil  section  is  provided.  Following 
that,  shown  in  Figure  28,  is  a  summary  of  the  major  geometric  characteristics 
of  such  planform;  note  that  fin  planform  geometry  data  is  given  for  one  panel 
of  each  fin  set,  since  it  is  assumed  that  each  fin  of  a  fin  set  is  identical.  If  a 
panel  is  made  up  of  multiple  segments,  the  geometric  data  is  provided  by 
panel  segment  (each  segment  is  assigned  a  number  starting  at  the  root).  Total 
panel  set  of  characteristics  is  also  provided.  This  total  panel  data  represents  an 
equivalent  straight-tapered  panel,  which  is  used  for  most  of  the  aerodynamic 
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calculations.  The  thickness-to-chord  ratio  shown  for  each  segment  is  that 
value  at  the  segment  root;  for  the  total  panel,  it  is  an  "effective"  value. 

If  an  airbreathing  inlet  is  specified  the  output  is  similar  to  that  in 
Figure  29.  This  output  reflects  the  user  input  definition  for  the  inlet  design 
specified.  It  is  provided  if  the  PRINT  GEOM  INLET  or  PART  control  cards  are 
included  in  the  input  case. 

4.2.2  Aerodynamic  Partial  Output 

The  output  on  the  configuration  aerodynaniics  is  most  extensive  when 
PRINT  AERO  or  PART  is  specified.  Output  is  created  for  the  body  and  each  fin 
set  on  the  configuration.  In  addition,  for  any  subsonic/ transonic  Mach 
number  (less  than  1.4)  an  analysis  by  the  Airfoil  Section  Module  is  made, 
which  involves  a  potential  low  analysis  of  the  airfoil  section  using  conformal 
mapping.  If  a  conhguration  has  inlets  additional  partial  output  is  included  to 
summarize  the  inlet  external  aerodynamics. 

If  base-jet  plume  interaction  calculations  are  specified  (BASE=.TRUE 
in  namelist  A^BOD),  then  there  will  be  one  or  two  separate  pages  of  output. 
Figure  30  shows  an  example  of  the  first  page  of  output.  This  page  will  always 
be  printed  if  BASE=.TRUE.  The  base  pressure  coefficient,  axial  force 
coefficient,  and  freestream  pressure  and  temperature  ratios  are  shown  versus 
angle  of  attack.  Also,  the  incremental  forces  and  moments  due  to  separation 
are  shown  versus  angle  of  attack.  If  extrapolation  of  the  base  pressures  and 
separation  conditions  database  occurs,  a  warning  message  is  printed  explaning 
what  input  variable  required  extrapolation.  A  second  page  of  output 
containing  the  boattail  separation  parameters  will  be  printed  if  there  are  any 
fins  on  the  missile  boattail.  The  separation  location  aft  of  the  nose  and  the 
Mach  cone  angle  are  shown  versus  angle  of  attack  for  each  panel  on  the  fin 
set.  Figure  31  shows  an  example  of  this  page.  This  output  is  provided  if  the 
PRINT  AERO  BODY  or  PART  control  card  is  input. 

The  protuberance  partial  output  is  printed  if  PRINT  AERO  BODY  or 
PART  is  used.  This  output  will  only  be  shown  if  the  namelist  PROTUB  is 
present  in  the  input  file.  Figure  32  is  an  example  of  the  protuberance  output. 
Protuberance  type,  location,  number,  and  axial  force  coefficient  are  listed  for 
each  protuberance  set.  The  total  axial  force  coefficient  or  zero  lift  drag 
coefficient  is  printed  at  the  bottom  of  the  page. 

As  shown  in  Figure  33,  the  body  alone  partial  aerodynantic  output  for 
normal  force  lists  the  axial  force  contributors,  potential  normal  force  (CN- 
POTENTIAL),  viscous  normal  forces  (CN-VISCOUS),  potential  pitching 
moment  (CM-POTENTIAL),  viscous  pitching  moment  (CM-VISCOUS),  and 
the  crossflow  drag  coefficient  (CDC).  The  cross-flow  drag  proportionality  factor 
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at  subsonic  and  transonic  speeds  is  also  given  for  reference.  These  data  are 
similar  to  that  obtained  for  elliptical  bodies. 

Figure  34  details  the  fin  normal  force  calculations  by  fin  set.  Each 
panel’s  contribution  to  the  configuration  normal  force  is  described.  The 
column  titled  CN-POTENTIAL  is  the  potential  contribution  and  the  column 
titled  CN-VISCOUS  is  the  viscous  contribution.  Their  sum  is  given  in  the 
column  titled  CN-TOTAL.  CNAA  is  the  nonlinear  variation  of  normal  force 
due  to  angle  of  attack  and  ALPHA  EQUIV  is  the  panel  angle  of  attack  due  to 
its  roll  position  on  the  body.  Figure  35  illustrates  the  fin  axial  force 
contributors  and  Figure  36  presents  an  example  of  the  fin  pitching  moment 
contributors. 

The  analysis  by  the  Airfoil  Section  Module  is  provided  in  a  format 
similar  to  Figure  37.  If  any  Mach  number  specified  produces  supersonic  flow 
on  the  airfoil  surface,  the  message  "CREST  CRITICAL  MACH  NUMBER 
EXCEEDED "  will  be  printed;  approximation  of  the  airfoil  section  data  is  then 
assumed.  These  fin  aerodynamic  increments  are  repeated  for  each  fin  set  on 
the  configuration.  Note  that  the  Airfoil  Section  Module  assumes  that  the 
panels  have  sharp  trailing  edges.  Any  panel  input  with  a  non-sharp  trailing 
edge  will  have  its  aerodynamic  characteristics  set  as  though  the  airfoil  was 
"ideal".  This  assumption  is  approximate  for  preliminary  design. 

Figure  38  shows  the  aerodynamic  output  available  when  inlets  are 
specified  on  the  configuration.  It  is  provided  when  PRINT  AERO  INLET  or 
PART  is  specified  in  the  case  inputs.  The  aerodynamics  summarized  for  inlets 
can  include  additive  drag  results  if  the  user  input  the  additive  drag 
calculation  flag.  The  maximum  mass  flow  ratio  is  printed  at  the  bottom  of  the 
page  if  the  additive  drag  is  calculated.  If  additive  drag  cannot  be  calculated,  a 
warning  message  is  printed. 

After  the  aerodynamic  details  for  each  component  of  the  configuration 
are  output,  the  aerodynamic  calculations  for  the  synthesis  of  the  complete 
configuration  follows.  For  the  example  case,  fin  set  1  results  would  be 
followed  by  fin  set  2  results  for  each  of  the  following  outputs: 

•  "FIN  SET  PRESENCE  OF  THE  BODY "  -  This  summarizes  the 
aerodynamic  incrementals  of  the  most  forward  set  of  fins 
with  the  influence  of  the  body.  Figure  39  presents  the 
example  of  this  output.  The  left-most  six  columns  include 
the  effect  of  body-on-fin  component  interference.  The  right¬ 
most  columns  represent  the  contribution  to  each  panel  to 
configuration  aerodynamics,  and  include  the  effect  of  body- 
on-fin  interference,  these  values  are,  in  effect,  individual 
panel  loads.  The  panel  characteristic  values  included  are 
"AEQn"  (the  panel  equivalent  (local)  angle  of  attack)  and 
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"CNn"  (the  panel  normal  force  coefficient).  The  sign 
convention  is  as  follov.s:  a  positive  panel  normal  force, 
hence,  equivalent  angle  of  attack,  produces  a  negative  roll 
moment.  Therefore,  panels  on  the  right  side  of  the 
configuration  will  produce  loads  and  angles  of  attack  opposite 
in  sign  to  those  on  the  left  side  of  the  configuration  even 
though  they  produce  the  same  physical  force  loading. 

•  "BODY-FIN  SET"  -  Aerodynamics  for  the  body  plus  most 
forward  set  of  fins  configuration.  It  is  produced  through 
addition  of  the  body  alone  and  wing  in  presence  of  the  body 
incrementals,  described  above.  The  results  shown  in  Figure 
40,  include  the  component  carryover  factors  K-W{B)  (wing  in 
presence  of  the  body  carryover  due  to  angle  of  attack),  K-B(W) 
(body  in  presence  of  the  wing  carryover  due  to  angle  of 
attack),  KK-B(W)  (body  in  presence  of  the  wing  carryover  due 
to  panel  deflection),  XCP-W(B)  (wing  in  presence  of  the  body 
carryover  center  of  pressure),  and  )(CP-B(W)  (body  in 
presence  of  the  wing  carryover  center  of  pressure).  This 
output  is  repeated  for  the  body  plus  each  additional  aft  fin  set, 
if  one  exists  on  the  configuration.  This  example  includes  two 
fin  sets  so  the  next  page  of  partial  output  would  look  like 
Figure  41.  If  additional  fin  sets  are  present  on  the 
configuration  additional  pages  are  output  with  each  one 
successively  included. 

•  "CARRYOVER  INTERFERENCE  FACTORS"  -  This  page  of 
partial  output  summarizes  the  carryover  factors  listed  in  the 
paragraph  above.  These  were  included  in  the  body  plus  fin  set 
calculations.  An  example  of  this  output  is  presented  in  Figure 
42. 


•  "COMPLETE  CONFIGURATION"  -  Complete  configuration 
aerodynamics.  This  output  is  illustrated  in  Figure  24.  The 
values  are  obtained  by  summing  the  body-wing  and  tail  in 
the  presence  of  the  wing  flow  field  data. 

In  addition  to  the  output  described  above,  more  data  is  presented  when 
the  BUILD  control  card  is  used.  Static  aerodynamics  are  output  for  each 
configuration  component.  Body  alone  aerodynamics  are  shown  in  Figure  43. 
Fin  alone  aerodynamics  are  shown  for  each  fin  set  present.  Figure  44  shows 
the  output  for  the  first  fin  set.  Static  aerodynamics  for  a  configuration  with 
body  plus  most  forward  set  of  fins  is  given  next.  Figure  45  shows  an  example 
of  this  output.  This  output  is  repeated  for  configurations  including  the  body 
plus  each  additional  fin  set  present. 
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If  the  PRINT  AERO  BEND  or  PART  control  card  is  used,  the  code  will 
compute  and  print  panel  bending  moment  coefficients  for  each  fin  set  on  a 
separate  page.  One  page  is  shown  in  Figure  46.  The  sign  convention  is  that 
assumed  for  the  individual  panel  loads  and  equivalent  angles  of  attack,  noted 
above.  The  bending  moment  coefficients  are  based  upon  the  reference  area 
and  longitudinal  length  given  at  the  top  of  the  page.  The  moments  are 
referenced  about  the  fin-body  structure  specified  by  the  root  chord  span 
station. 

Figure  47  illustrates  the  panel  hinge  moments  coefficients  computed 
when  the  control  cards  PRINT  AERO  HINGE  or  PART  are  used.  The 
reference  area  and  longitudinal  reference  length  given  at  the  top  of  the  page 
are  used-  All  moments  are  computed  about  the  hinge  line,  which  is  defined 
using  namelist  DEFECT. 

If  TRIM  is  specified,  the  user  can  selectively  print  the  six  untrimmed 
static  aerodynamic  tables  used  in  the  trim  process.  An  example  is  shown  in 
Figure  48.  The  code  computes  the  six-component  aerodynamics  at  ten 
deflection  angles  for  each  specified  angle  of  attack,  then  interpolates  for  Cm=0. 
Note  th’t  this  trim  process  can  be  used  to  create  contiol  authority  data, 
effectively  giving  the  user  10  deflection  angles,  20  angles  of  attack,  and  20 
Mach  numbers  per  input  case. 

4.2.3  Pressure  Distribution  Data 

If  the  Mach  number  is  supersonic  (M  S  1.2),  the  user  hab  the  option  to 
print  the  surface  pressure  distributions  over  the  body  and  fins.  This  option  is 
selected  only  through  the  addition  of  the  control  card  PRESSURES.  Since 
three  body  alone  supersonic  methods  are  available  (Van  Dyke  Hybrid, 
Second-Order  Shock  Expansion  (SOSE),  and  Newtonian  flow)  the  capability 
exists  to  output  the  pressure  distribution  data  from  any  one  of  these  methods. 
The  method  to  be  used  in  the  calculation  of  the  press.  ' e  data  is  controlled 
with  the  control  cards  SOSE  and  HYPER;  if  neither  consul  card  is  input,  the 
Van  Dyke  Hybrid  method  is  selected.  Because  of  the  nature  of  the 
calculations,  body  alone  pressures  are  printed  for  angles  of  attack  less  than  or 
equal  to  15  degrees  when  using  the  Hybrid  or  SOSE  techniques. 

The  capability  also  exists  for  the  user  to  output  the  pressure 
distribution  data  over  fins  at  any  Mach  number  greater  than  1.05.  This  option 
is  also  controlled  by  the  PRESSURES  control  card.  Due  to  the  nature  of  the 
method,  only  pressure  distribution  data  at  zero  angle  of  attack  is  presently 
output. 

Figures  49,  50,  and  51  illustrate  typical  output  produced  when 
PRESSURES  is  specified.  The  format  of  Figure  49  is  only  available  when  SC)SE 
is  specified;  all  other  body  alone  pressure  methods  produce  output  similar  to 
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Figure  50  for  bodies.  Figure  51  is  representative  of  fin  pressure  distribution 
output.  Note  that  calculation  of  pressures  is  a  time-consuming  process;  much 
higher  computational  times  will  be  required. 

All  body  pressure  distribution  data  is  based  on  a  configuration  that  has 
body  diameter  of  unity;  that  is,  the  configuration  is  expressed  in  calibers  (or 
body  diameters).  The  longitudinal  stations  at  which  pressure  coefficient  data 
is  desired  cannot  be  usrj  specified;  however,  sufficient  data  is  provided  to 
permit  accurate  interpolation  for  most  applications. 

4.3  Dynamic  Deriva  ''es 

As  shown  in  Figure  52,  the  total  configuration  dynamic  derivatives  are 
provided  in  compact  form  for  easy  interpretation.  The  dynamic  derivatives 
are  summarized  as  a  function  of  angle  of  attack  in  the  user  specified  units. 
The  coefficients  provided  are  as  follows: 

CNQ  Normal  force  coefficient  due  to  pitch  rate 

CNAD  Normal  force  coefficient  due  to  acceleration  in 
angle  of  attack  (a) 

CMQ  Pitching  moment  coefficient  due  to  pitch  rate 

CMAD  Pitching  moment  coefficient  due  to  acceleration  in 
angle  of  attack  (a) 

Note:  For  body  alone  and  body  +  fin  set  data  CMQ  and  CMAD  are 
pre«eided  as  the  sum  CMQ+CMAD. 

The  dynamic  derix  atives  are  printed  after  all  static  coefficients  and 
partial  static  aerodynamics  are  printed.  If  a  BUILD  or  PART  card  is  input, 
additional  dynamic  derivatives  for  partial  configurations  and/or 
configuration  components  are  printed. 

4  External  Data  Files 


The  code  has  the  capability  to  be  used  in  conjunction  with  other 
missile  design  tools,  such  as  post-processing  plotting  programs  or  trajectory 
programs.  Fixed  format  aerodynamic  data  is  output  as  an  external  data  file 
with  the  addition  of  the  PLOT  control  card.  Included  in  this  data  file  are  the 
six  component  forces  and  moments  based  upon  the  user  specified  reference 
quantities.  In  order  to  print  component  buildup  data  to  the  plot  file  the 
BLULD  and  PLOT  control  cards  must  be  present  in  the  case. 
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An  option  to  create  a  user  specified  format  data  file  is  a.  jo  available. 
The  control  cards  WRITE  and  FORMAT  have  been  designed  for  easy  access  to 
this  capability. 

4.5  Extrapolation  Messages  and  Array  IXimps 

As  shown  in  Figure  53,  the  extrapolation  messages  are  summarized  for 
all  design  charts  which  have  been  extrapolated  during  the  execution  of  the 
case.  Since  many  of  the  aerodynamic  methods  do  not  include  design  charts, 
but  are  either  closed-form  equations  or  complete  theoretical  methods,  this 
option  is  most  useful  in  the  subsonic  and  transonic  Mach  regimes. 
Extrapolation  messages  are  only  provided  if  the  control  card  PRINT  EXTRAP 
appears  in  the  case  inputs.  The  data  titled  "ROUTINE  TRACE-BACK"  lists  the 
subroutines  called  when  the  look-up  was  performed;  "X"  represents  the 
independent  variable  and  "Y"  represents  the  dependent  variable  in  the 
extrapolation. 

When  it  is  necessary  to  examine  the  values  stored  ir  internal  data 
arrays  the  DUMP  control  card  can  be  used.  This  control  card  causes  the 
contents  of  the  named  data  arrays  to  bo  printed  in  a  form  similar  to  Figure  54. 
Array  dumps  are  provided  for  each  Mach  number  of  the  input  case,  and 
represent  the  data  block  contents  at  aerodynamic  calculation  completion. 

Note  that  all  data  arrays  are  initialized  to  a  constant  named 
"UNUSED",  which  is  preset  to  a  value  of  1  x  Hence,  any  array  element 
which  contains  this  constant  was  not  changed  during  execution  of  the  case 
(since  it  is  highly  unlikely  that  this  constant  will  result  from  any  calculation). 
This  scheme  permits  rapid  "tracking"  of  program  calculation  sequences  while 
in  "debug"  mode. 
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Figure  21  Input  Error  Checking  Output 
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Figure  22 


Case  Input  Listing 
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Figure  23  Example  of  Default  Substitutions 
for  Incomplete  Case  Inputs 
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Figure  24  Total  Configuration  Output  Summary 


78 


LIRAtt  CATA  rOR  BCOY  AI/>R  «A8  ClRAArCD  USIWC  ?BS  8tCCaX>-0ia>Ul  SHOCK  IX?AH8ZOV  MCTaOD 


E 


c 


t; 


3 


g  S 

&  §  S  . 

o  o  |l 


<>»  »>l 


Figure  25  Trimmed  Output  Summary 
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Figure  26  Body  Ge<)meiry  Output 
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Figure  28  Fin  Geometry  Output 
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Figure  29  Inlet  Geometry  Output 
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Figure  30  Base-Jet  Plume  Interaction  Output  -  Page  1 
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Figure  31  Base-Jet  Plume  Interaction  Output  -  Page 
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Figure  32  Protuberance  Output 
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I'igure  33  body  Alone  Aerodynamic  Partial  Output 
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Figure  34  Fin  Normal  Force  Partial  Output 
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Figure  36  Fin  Pitching  Moment  Partial  Output 
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Figure  37  Airfoil  Section  Aerodynamic  Partial  Output 
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Figure  38  Inlet  Aerodynamic  Partial  Output 
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Figure  40  Body  Plus  Fin  Sel  Partial  Output 
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LIHLAR  CATA  FOR  BODY  ALOKX  WAS  GCWtRATCD  USIHC  THE  SBC CVD-ORDER  SKCXZR  EXPANSION  METHOO 
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Figure  45  Body  Pius  Fin  Static  Aerodynamic  Partial  Output 
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Figure  47  Panel  liinge  Moment  Partial  Output 
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figure  49  Body  Pressure  Distribution  from  SOSE,  AOA=0° 
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Figure  50  Body  Pressure  Distribution  at  Angle  of  Attack 
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Figure  51  Fin  Pressure  Distribution  Output 
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Figure  52  Dynamic  Derivative  Output 
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Figure  53  Extrapolation  Message  Output 
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Figure  54  Internal  Array  Dump  Output 


A.  EXA>vlPLE  PROBLEMS 


This  appendix  presents  two  example  problems  for  use  in  verifying  the 
proper  operation  of  the  computer  code,  or  for  use  as  a  model  for  the  setup  of 
inputs  of  similar  configurations.  The  first  example  is  a  simple  tangent  ogive 
nose-cylinder  circular  l^y.  It  has  a  planar  wing  (two  panels)  and  a  cruciform 
set  of  tails  orientated  in  the  "plus"  configuration.  There  are  two  input  cases 
for  this  problem.  The  first  demonstrates  the  output  generated  when  the 
PART  control  card  is  included.  The  second  demonstrates  the  output  created 
when  trim  has  been  requested;  note  that  trim  partial  output  has  been 
requested  with  the  PRINT  AERO  TRIM  control  card. 

The  second  problem  is  a  body-tail-inlet  configuration.  Three  Mach 
numbers  are  requested,  one  subsonic,  one  transonic,  and  one  supersonic. 
Although  all  three  Mach  numbers  could  be  run  in  a  single  case,  they  have 
been  divided  into  three  separate  cases  to  illustrate  the  "SAVE"  feature  as  well 
as  selecting  particular  output  for  illustration. 
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A.l  Example  Problem  1 


The  first  example  problem  is  shown  in  Figure  A-1.  It  is  comprised  of  a 
3-caliber  tangent  ogive  nose  attached  to  a  cylindrical  body;  a  triangular 
monoplane  set  of  wings;  and  a  cruciform  set  of  tails  orientated  in  the  "plus" 
position.  The  first  case  is  a  simple  angle  of  attack  sweep;  component  buildup 
data  and  partial  output  are  requested.  The  second  case  is  a  trim  of  the 
configuration  using  the  two  horizontal  tail  surfaces.  The  inputs  are  shown  in 
Figure  A-2. 
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Figure  A-2  Example  I’roblem  1  Input 
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A. 2  Example  Problem  2 


The  configuration  for  this  example  is  sketched  in  Figures  A-3  and  A-4. 
The  figure  is  a  modified  copy  of  the  wind  tunnel  model  drawings  from 
NASA  Technical  Memorandum  84557.  The  model  definition  in  these  figures 
is  representative  of  the  detail  normally  found  on  design  drawings. 

This  example  has  been  divided  into  subsonic,  transonic,  and 
supersonic  cases.  Each  case  is  run  for  one  Mach  number.  Although  all  three 
Mach  numbers  could  have  been  run  in  one  case  they  were  run  separately  to 
demonstrate  the  SAVE  capability. 

This  example  provides  a  check  case  for  the  inlet  option.  It  can  be  used 
by  the  user  to  make  sure  that  he  understands  the  inputs.  Figure  A-5  shows 
the  inputs  required  to  run  this  example. 
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Figure  A-3  Example  Problem  2  Configuration 


INLET-INLET  FAIRING  DETAILS 
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Figure  A-4  Example  Problem  2  Tail /Inlet  Details 


Figure  A-5  Example  Problem  2  Input 
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Figure  A-5  Example  Problem  2  Input  (Continued) 
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B.  PLOT  FILE  FORMAT 


When  the  PLOT  control  card  is  used,  a  formatted  data  file  is  written  to 
unit  3  which  can  be  used  in  a  separate  plotting  program.  The  file  produced 
has  the  following  format: 


1  LINE 

COLUMN 

CONTENT 

1 

K3 

Word  RUN 

4-7 

A  sequential  run  number  beginning  with 
one  (1)  and  incrementing  by  one  (1) 

8-10 

Number  of  angles  of  attack 

13 

Units  System: 
l-inches 

F-fect 

C-centimeters 

M-meters 

14-23 

I.D.  Field  (sec  below) 

2 

1-10 

Mach  Number 

11-20 

Reynolds  Number/foot  (or  meter) 

3 

1-10 

Reference  Area 

11-20 

Longitudinal  Reference  Length 

21-30 

Longitudinal  C.C.  Location  (X  c^) 

31-40 

Lateral  Reference  Length 

41-50 

Vertical  C.G.  Location  (Zcg) 

4 

(Repeated  for 
each  angle  of 
attack) 

1-10 

Angle  ot  Attack 

11-20 

Cn  or  ajRiM 

21-30 

Cm  or  CntRIM 

31-40 

Ca  orCATRIM 

41-50 

Cy  u.  CYi'KlM 

1-60 

Cn  OT  CnTRlM 

1-70 

ClorCiTRIM 

LAST 

1 

Character  R  (denotes  end  of  data  for  set) 

The  I.D.  field  identifies  the  configuration  for  which  the  data  is 
applicable. 


ID.  Code 

Configuration/Condition 

B 

Body  Alone 

FI 

Most  forward  fin  set  or  fin  set  #1 

F2 

Second  most  forward  fin  set  or  fin  set  #2 

F3 

Third  most  forward  fin  set  or  fin  set  #3 

F4 

Fourth  most  forward  fin  set  or  fin  set  #4 

BFl 

Body  plus  most  forward  fin  set 

BF12 

Body  plus  two  most  forward  fin  sets 

BFl  3 

Body  plus  three  most  forward  fin  sets 

BF14 

Body  plus  four  most  forward  fin  sets 

Dn 

Untrimmed  data  for  deflection  angle  #n  (n=l  to  10) 

TRIMMED 

Trimmed  results 
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Component  buildup  data  is  written  to  the  file  if  the  BUILD  control  card 
is  used.  If  the  TRIM  control  card  is  used,  both  the  trimmed  and  untrimmed 
results  are  output.  Otherwise,  only  the  final  configuration  data  output. 


